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PRÓLOGO 

Este  trabajo  ha  surgido  como  consecuencia  de  explicar  durante  varios  años 
el  tema  relativo  al  estudio  de  los  métodos  numéricos  para  la  resolución  de  ecua¬ 
ciones  escalares  no  lineales,  así  como  de  una  primera  idea  que  apareció  en  el 
desarrollo  de  dicho  tema,  la  cual  consiste  en  relacionar  la  variación  de  la  conve¬ 
xidad  de  la  curva  y  =  f(x)  con  la  velocidad  de  convergencia  que  un  proceso  itera¬ 
tivo  tiene  al  aproximar  una  raíz  de  la  ecuación  f(x)=0. 

A  partir  de  esta  primera  idea,  y  junto  a  la  Dra.  M.  A.  Salanova,  comenzamos 
a  realizar  algunos  trabajos  de  investigación  sobre  las  ecuaciones  escalares  no 
lineales  como  paso  previo  al  estudio  de  la  resolución  de  operadores  no  lineales 
en  espacios  de  Banach,  que  era  nuestro  principal  objetivo.  Es  como  resultado  de 
esta  investigación  previa  que  hemos  elaborado  estas  notas. 

Finalmente  queremos  expresar  nuestro  agradecimiento  a  nuestros  compa¬ 
ñeros  de  la  Universidad  de  La  Rioja  que  a  lo  largo  de  estos  años  nos  prestaron 
su  inestimable  ayuda.  En  particular,  queremos  agradecer  a  los  profesores  J.  A. 
Ezquerro  y  J.  M.  Gutiérrez  sus  observaciones  e  indicaciones  en  la  correción  y 
mejora  de  estas  notas. 


Logroño,  Febrero  de  1996. 


Miguel  Angel. 


ÍNDICE 

TEMA  I.  ÍNDICES  DE  MEDIDA  DE  LA  CONCAVIDAD  Y 
CONVEXIDAD  DE  UNA  CURVA 

INTRODUCCIÓN .  13 

CAPÍTULO  I.  Medida  de  la  convexidad  de  una  curva 
mediante  el  operador  logaritmo 

1 .0.  Introducción  .  15 

1.1.  Funciones  logarítmicamente  convexas  .  19 

1 .2.  Definición  de  grado  del  convexidad  logarítmico  de  una  función  ....  23 

1 .3.  Propiedades  del  grado  de  convexidad  logarítmico .  34 

CAPÍTULO  II.  Medida  de  la  concavidad  de  una  curva  mediante 
el  operador  exponencial 

2.0.  Introducción  .  37 

2.1.  Funciones  exponencialmente  cóncavas  .  38 

2.2.  Definición  del  grado  de  concavidad  exponencial  de  una  función  ...  41 

2.3.  Propiedades  del  grado  de  concavidad  exponencial  .  47 

CAPÍTULO  III.  Indices  de  medida  de  la  convexidad  y 
de  la  concavidad  de  una  función 

3.0  Introducción  .  49 

3.1.  Funciones  cp  -  convexas  .  51 

3.2.  rp  -  grado  de  convexidad  .  54 

3.3.  Grado  de  convexidad  .  62 

3.4.  Grado  de  concavidad  .  63 


9 


TEMA  II.  LA  CONCAVIDAD  Y  CONVEXIDAD  EN  LA  CONVERGENCIA 
DE  PROCESOS  ITERATIVOS 


INTRODUCCIÓN  .  67 


CAPÍTULO  I  La  concavidad  y  la  convexidad  en  la  convergencia 
de  procesos  iterativos  de  orden  menor  que  tres  . 

1.1.  El  método  de  Whittaker .  69 

1.1.1.  Generalidades  .  69 

1 .1 .2.  La  concavidad  en  el  método  de  Whittaker .  71 

1 .2.  El  método  de  Newton  -  Raphson  .  78 

1.2.1.  Generalidades .  78 

1 .2.2.  La  convexidad  en  el  método  de  Newton-Raphson .  80 


CAPÍTULO  II.  La  concavidad  y  la  convexidad  en  la  convergencia 

de  procesos  iterativos  para  la  resolución  de  ecuaciones 
no  lineales 

2.0.  Introducción  .  87 

2.1.  El  método  de  Chebyshev  .  87 

2.1 .1.  Generalidades  .  87 

2.1 .2.  Discusión  de  la  convergencia .  89 

2.1.3.  Notas  prácticas .  92 

2.2.  El  método  de  Halley  .  95 

2.2.1.  Generalidades .  95 

2.2.2.  Discusión  de  la  convergencia .  96 

2.2.3.  Notas  prácticas .  98 

2.2.4.  Apéndice . 101 

2.3.  Familia  de  procesos  iterativos  de  tercer  orden  . 103 

2.3.1.  Discusión  de  la  convergencia . 103 

2.3.2.  Notas  prácticas . 106 

BIBLIOGRAFÍA  . 109 


10 


INTRODUCCIÓN 

El  estudio  de  la  concavidad  y  convexidad  de  un  curva  es  un  antiguo  proble¬ 
ma  resuelto  por  los  matemáticos.  Está  perfectamente  caracterizado  el  hecho  de 
que  una  curva  sea  cóncava  o  convexa.  Sin  embargo  el  problema  de  medir  esta 
concavidad  o  convexidad  es  una  cuestión  menos  desarrollada.  Son  conocidos  los 
grados  de  convexidad  de  una  curva  obtenidos  por  JENSEN  y  POPOVICIU  [3];  no 
obstante  dichos  grados  de  convexidad  tienen  un  eminente  sentido  teórico.  Su 
aplicación  práctica  resulta  excesivamente  complicada  pues,  entre  otras  cosas, 
conlleva  la  consideración  de  particiones  de  un  intervalo  según  diferentes  medi¬ 
das,  ya  que  estos  grados  de  convexidad  son  globales  sobre  dicho  intervalo. 

Otra  forma  de  medir  la  convexidad  de  una  curva  se  sugiere  en  el  Teorema 
de  BOHR-MOLLERUP  [2],  para  la  definición  de  la  función  Gamma.  En  este 
resultado  aparece  el  concepto  de  función  logarítmicamente  convexa  [19],  es 
decir,  una  función  cuyo  logaritmo  es  una  función  convexa.  J.  GARAY  y  M.  A. 
HERNANDEZ  [6],  introducen  el  grado  de  convexidad  logarítmico.  Mediante  la 
sucesiva  aplicación  del  operador  logaritmo,  operador  que  “concaviza”  a  una  fun¬ 
ción  convexa,  se  obtiene  una  medida  puntual  de  la  convexidad,  dada  por  la 
resistencia  de  esta  función  a  dicho  operador.  Es  decir,  el  número  de  veces  que 
hay  que  aplicar  el  operador  logaritmo  a  la  función  para  obtener  como  resultado 
una  función  cóncava  nos  da  una  medida  puntual  de  la  convexidad  de  la  función 
que  nos  permite  definir  el  grado  de  convexidad  logarítmico,  cuyo  estudio  se 
desarrolla  en  el  primer  tema. 

A  continuación,  en  el  segundo  tema,  apoyándonos  en  la  interpretación  geo¬ 
métrica  del  grado  de  convexidad  logarítmico  definimos  el  correspondiente  grado 
de  concavidad  exponencial  para  una  función  cóncava. 

Terminamos  nuestro  estudio  de  la  concavidad  y  convexidad  de  una  curva 
extendiendo  estos  conceptos  para  operadores  9  cualesquiera,  definimos  así  el 
9  -grado  de  convexidad  y  el  9  -grado  de  concavidad  de  una  función  convexa  y 
cóncava  respectivamente. 
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En  este  estudio  consideraremos  que  una  función  f:  (a,b)  R  se  denomina 
convexa  si  se  verifica  que  f(A  x  +  (1-A,)  y  )  <  X  f(x)  +  (1-A.)  f(y)  para  todo  x,y  e 
(a,b)  y  X  e(0,1).  Se  dirá  estrictamente  convexa  si  la  desigualdad  anterior  es  estric¬ 
ta  para  todo  x  *  y  en  (a,b).  Ahora  bien,  dado  que  en  nuestro  estudio  considera¬ 
remos  en  todo  momento  funciones  f  eC(P((a,b)),  p  >  2,  es  conocido  [19]  que  en 
esta  situación  f  es  convexa  sí  y  sólo  sí  f”(x)  >  0.  Y  si  f"(x)  >  0  en  (a,b)  enton¬ 
ces  se  verifica  el  mismo  resultado  para  f  estrictamente  convexa.  A  partir  de  esta 
caracterización,  dado  que  en  nuestro  estudio  intentamos  definir  una  medida  pun¬ 
tual  de  la  convexidad  de  una  función  consideraremos  la  convexidad  de  forma 
puntual.  Dado  el  carácter  local  que  tiene  la  convexidad,  es  difícil  establecer  una 
definición  de  función  convexa  en  un  punto.  Teniendo  en  cuenta  la  restricción  que 
ello  supone,  nosotros  consideraremos,  para  nuestro  estudio,  puntos  x0  donde  la 
función  sea  convexa  verificando  que  f”(xo)>0,  en  cuyo  caso  diremos  que  f  es  de 
carácter  estrictamente  convexo  en  x0,  o  bien  puntos  x0  tales  que  existe  ke  N, 
número  par,  con  k  <  p,  tal  que  f”(x0)  =...  =  f^'^(x0)  =  0  y  f^(x0)  >  0.  En  esta 
situación  diremos  que  f  es  de  carácter  convexo  en  sentido  no  estricto  en  x0. 
Cualquier  otra  situación  de  convexidad  en  un  punto  x0  queda  fuera  de  nuestro 
estudio.  Notemos  que  la  función  x^  es  convexa  en  sentido  estricto  en  el  punto 
xo=0,  sin  embargo  según  nuestra  consideración  es  de  carácter  no  estricto  en 
dicho  punto.  En  el  caso  de  la  concavidad  consideramos  la  situación  análoga.  Este 
mismo  criterio  lo  seguiremos  con  los  puntos  críticos. 
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CAPÍTULO  I 

MEDIDA  DE  LA  CONVEXIDAD  DE  UNA  CURVA  MEDIANTE 
EL  OPERADOR  LOGARITMO 


1.0.  Introducción 

El  Teorema  de  BOHR-MOLLERUP,  relativo  a  la  definición  de  la  función 
gamma  [2],  utiliza  el  hecho  de  que  dicha  función  es  “logarítmicamente  convexa” 
en  R+,  es  decir,  su  logaritmo  es  una  función  convexa,  (fig.1.) 
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Resulta  sencillo  observar  que  esta  propiedad  no  es  verificada  por  todas  las 
funciones  convexas.  Basta  considerar,  por  ejemplo,  la  función  x^.  (fig.2.) 


Esta  situación  nos  lleva  a  pensar  que,  en  cierto  sentido,  existen  funciones 
“más  convexas”  que  otras.  De  hecho,  lo  que  estamos  considerando  es  la  acción 
del  operador  logaritmo  sobre  el  conjunto  de  las  funciones  convexas.  Acción  que, 
por  ser  el  logaritmo  una  función  cóncava,  “concaviza”  las  funciones  a  las  que  se 
les  aplica.  En  este  capítulo  vamos  a  obtener  un  índice  de  medida  puntual  de  la 
convexidad  de  una  función.  Para  lograr  esto,  mediremos  la  resistencia  que  ofre¬ 
ce  una  función  convexa  a  dejar  de  serlo  mediante  sucesivas  aplicaciones  del  ope¬ 
rador  logaritmo,  es  decir,  calcularemos  el  número  de  veces  que  debemos  aplicar 
dicho  operador  a  la  función  convexa  considerada,  para  que  el  resultado  deje  de 
ser  una  función  convexa. 
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La  figura  3  muestra  que  hemos  de  aplicar  tres  veces  el  operador  logaritmo 
para  que  la  función  f(x)  =  x  -  eos  x  +  4,  que  es  convexa  en  x0  =  0,  deje  de  serlo. 

Este  número  nos  proporcionará  un  índice  de  medida,  en  cada  punto,  de  la 
convexidad  de  la  función  considerada. 

Notemos  por  otra  parte  que  si  consideramos  una  función  f  en  las  condiciones 
indicadas,  normalizada  por  f(x0)  =  1  para  simplificar  nuestros  cálculos,  entonces 
tenemos  que  la  curvatura  de  f  es  mayor  que  la  del  log  f  (fig.4.),  ya  que 

|r[iogf(x)]x=Xo  =  f'(x0)  y  ¿iflogí(x)]x=Xo  =  fw(x0)-f'(x0)2I 
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y  las  curvaturas  de  f(x)  y  de  log  f(x)  en  el  punto  x0  vienen  dadas  por 


[K(f(x))]x 


f"(x0) 


°  (1  +f'(x0))3/2 


[K(log  f  (x))] 

x  xo 


f"(x0)-  Hx/ 

3/  2 

(1  +f'(x0)) 


luego  parece  claro  que  se  verifica  nuestra  idea  geométrica  de  construcción  de 
éste  índice  de  convexidad. 
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1.1.  Funciones  logarítmicamente  convexas. 

Comenzamos  nuestro  estudio  analizando  el  concepto  de  función  logarítmica¬ 
mente  convexa  introducido  por  BOHR  -MOLLERUP.  Este  análisis  tiene  como  obje¬ 
tivo  el  de  conseguir  una  expresión  analítica  que  caracterice,  de  alguna  manera, 
dicho  concepto. 

Sea  V  un  entorno  del  punto  x0  y  f  una  función  derivable,  al  menos  de  orden 
dos  (f  e  c(m(V)  con  m  >  2  ),  positiva  y  convexa  en  V.  Sabemos  que  f  es  una  fun¬ 
ción  logarítmicamente  convexa  en  x0  si  log  f  es  una  función  convexa  en  x0. 
Entonces,  si  f  es  una  función  logarítmicamente  convexa  en  x0,  existe  un  entorno 
W  de  x0  y  una  función  h  e  c(m  (W)  convexa,  W  c  V,  de  manera  que  f(x)  =  exp 
[h(x)j  en  W,  es  decir,  h(x)  =  log  f(x)  en  W.  Operando  con  la  expresión  obtenida 
para  la  función  f  se  sigue  que 

f’(x)  =  h’(x)  exp  [h(x)j  =  h’(x)  f(x) 

f”(x)  =  h”(x)  f(x)  +  h’(x)  f’(x)  =  [h”(x)  +  h’(x)2]  f(x) 

y  de  estas  igualdades  se  obtiene 

(1.1)  h"(x)  =  lí-xlf-Tx)  -  [ f'(x)]2 

[f  (X)] 

igualdad  que  nos  va  a  permitir  realizar  el  estudio  de  la  convexidad  de  la  función 
h  y  por  tanto  de  log  f,  al  menos  en  el  caso  de  que  esta  convexidad  sea  en  senti¬ 
do  estricto. 

A  partir  de  (1.1),  considerando  su  valor  en  el  punto  x0  y  los  diferentes  tipos 
de  convexidad  que  se  pueden  presentar  en  él,  obtenemos  los  siguientes  resulta¬ 
dos,  diferenciando  el  caso  de  que  x0  sea  un  punto  crítico. 


Teorema  1.1. 

Sea  x0  un  punto  no  crítico  de  f 

(i)  Si  f  es  una  función  convexa  en  sentido  estricto  se  verifica 

(a)  f  es  logarítmicamente  convexa  en  x0  en  sentido  estricto,  es  decir,  log 
f  es  estrictamente  convexa  en  xQ  sí  y  sólo  sí 

[f'íxj]2 
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(b)  f  es  logarítmicamente  convexa  en  xQ  en  sentido  no  estricto,  sí  y  sólo 
sí  existe  keN,  número  par,  k<m,  de  manera  que 


(1.2) 


[f  (x  0)]1~Vt(x0) 
[f'ÍXo)]’ 


2  <t  <  k-  1 


,(k.  W  /  vk-1 
f  (Xp)f  (Xq) 

[f  (x^]" 


>  1 


(ii)  Si  f  es  una  función  convexa  en  sentido  no  estricto  entonces  f  no  es  loga¬ 
rítmicamente  convexa  en  x0. 

Demostración: 

En  cuanto  al  apartado  (i),  por  (1 .1)  el  punto  (a)  resulta  obvio.  Por  otra  parte, 
para  probar  (b),  como  h  =log  f  es  convexa  en  x0,  pero  en  sentido  no  estricto,  se 
verifica  que  existe  keN,  número  par,  de  manera  que  h”(x0)  =...  =  h(k‘1(x0)  =  0  y 
h(k(x0)  >  0,  luego  teniendo  en  cuenta  (1.1)  se  obtiene  que  f(x0)f”(x0)=  [f’(x0)]2. 

Para  probar  (1 .2)  consideramos  que  f’(x0)  =  h’(x0)  f(x0)  y  por  tanto  se  verifica 


(1.3)  f  (n+1(x ^  ^ ij")  h(t  +1(Xo)f(n  *  (x q) 

entonces 

n  =  0  f  (xQ)  =  h’(x0)  f(x0) 


n  =  1 


f"(x0)  =  h'(x0)f'(xo)  +h"(x0)f(xo) 


f'(x0) 

,(Xo>=W 


n  =  2  f'"(xo)  =  h'(xo)f"(xo)  +  L  Jh"(x0)r(x0)  +h"'(x0)f(xo) 


f... 


(xo) 


Hx/ 

f(xo)2 


k-2 


(k-2  f'(Xo) 

Suponiendo  que  f  (xJ  = - ¡— t- 

í  (x  o)  _ 


por  recurrencia  se  tiene 


.  „  ,  (k— 1  .  k-2fk-2V  (t+1.  (k-2-t  ,  ...  (k-2 

n  =  k  -  2  fv  (x1=  X  I  hv  (x  g)f '  (x0)  =  h(xo)f  (*o> 

t=0v  J 


f(k  ’(x0)  = 


f(x0) 


.  k-1 


f(xo) 


k-2 
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Por  otra  parte,  si  consideramos 

.  .  f(k  tVk  - 1\  (t+1  (k— 1— t  w(k-1  \  JK 

n  =  k -  1  f  (xQ)=  X^t  Jh  (x^f  (xQ)  =  h(x0)f  (xj  +  h  (xQ 

=>  f(k<xo)=  f  <X¿  +h(l<(xo)f(xo)  =» 
f(Xo) 

=>  f(k(x0)  fíXo)^1  -  f’(x0)k  =  h(k(x0)  f(x0)k. 

Como  k  es  par,  f(x0)  >  0  y  h(k(x0)  >  0  entonces  se  sigue  que  f(k(x0)  fíXQjk'1  > 

f’(x0)k- 

Para  ver  que  las  condiciones  también  son  suficientes,  basta  despejar  hí^Xg) 
en  cada  una  de  las  igualdades  anteriores  y  se  comprueba  fácilmente  que  h  es 
convexa  en  sentido  no  estricto. 


Para  el  apartado  (ii),  se  deduce  de  la  expresión  (1.1)  que  h”(x0)<  0.  Por  tanto 
h(x)  =  log  f  no  es  convexa  y  por  tanto  f  no  es  logarítmicamente  convexa  # 
Notemos  que  la  condición  (b)  del  apartado  (i)  para  t  =  2,  no  es  suficiente 
para  que  f  sea  logarítmicamente  convexa  en  x0  en  sentido  no  estricto,  ya  que  si 
f(x0)  f”(x0)[f’(x0)]  2  =  1 ,  no  podemos  asegurar  que  la  función  h=log  f  sea  convexa 
(fig.5.) 


21 


LA  CONVEXIDAD  EN  LA  RESOLUCIÓN  DE  ECUACIONES  ESCALARES  NO  LINEALES 


A  continuación  indicamos  unos  ejemplos  que  corroboran  este  hecho. 

Ejemplo  1.1. 

Sean  las  funciones 
f-l  (x)  =exp  [x  +  x4/4  ] 

Í2(x)  =exp  [x  -  x4/4  ] 

Í3(x)  =exp  [x  +  x3/3  ]  ( fig.5.  ) 

verifican  que  f¡(0)f”¡(0)  =  f’¡(0)2  i  =  1,2,3,  sin  embargo  h-|=logf-|  es  convexa  en  el 
origen,  h2=log  f2  es  cóncava  y  h3=log  Í3  tiene  un  punto  de  inflexión.  Esto  prueba 
la  indeterminación  que  presenta  esta  situación. 

Teorema  1 .2. 

Si  xQ  es  un  punto  crítico  de  f,  es  decir  en  nuestras  condiciones  un  mínimo, 
entonces  se  verifica 

(1.4)  lím.  =  +  ~ 

x^xo  [f  (X)] 

Además  x0  será  un  mínimo  de  h  =log  f,  manteniéndose  para  h  el  carácter 
estricto  o  no  estricto  que  x0  tenga  para  f. 

Demostración: 

Para  probarlo  es  necesario  distinguir  dos  situaciones  posibles.  A  parte  de 
verificarse  que  f’(x0)  =  0,  puede  ocurrir  que 

1 )  f”(x0)  >  0 

2)  Exista  keN,  número  par,  k  <  m,  de  manera  que 
f”(x0)  =■■■  =  f(k_1  (x0)  =  0  y  f(k(xQ)  >  0. 

En  el  caso  1 .)  resulta  obvio  que,  como  f(x0)  >  0  y  f”(x0)  >0,  se  verifica  (1 .4). 
Además,  teniendo  en  cuenta  (1.1),  se  tiene  que  h”(x0)  >  0,  siendo  h’(x0)  =  0, 
luego  se  mantiene  el  carácter  estricto  del  mínimo  x0  de  f  para  la  función  h  =log  f. 

En  cuanto  a  2.),  aplicando  la  regla  de  L'Hópital  [20],  se  prueba  que  también 
verifica  (1 .4).  Y  considerando  la  expresión  (1 .3)  se  obtiene  que  h’(x0)  =  h”(x0)  =... 
=  h(k_1  (x0)  =  0  y  h(k(x0)>  0.  # 

Pueden  estudiarse  algunas  propiedades  interesantes  de  las  funciones  loga¬ 
rítmicamente  convexas  en  el  libro  de  Roberts  [19]. 
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1.2.  Definición  del  grado  de  convexidad  logarítmico  de  una  función 

En  este  párrafo  vamos  a  establecer  un  índice  de  convexidad  de  una  función. 
Para  ello,  como  ya  indicamos  en  la  introducción  de  este  tema,  vamos  a  considerar 
el  operador  logaritmo.  Entonces,  reiterando  el  procedimento  utilizado  al  considerar 
funciones  logarítmicamente  convexas  visto  en  el  párrafo  anterior,  es  decir  la  aplica¬ 
ción  del  operador  logaritmo,  vamos  a  obtener  un  número  que  nos  va  a  permitir  esta¬ 
blecer  ese  índice  puntual  de  la  convexidad  que  indicábamos  anteriormente. 

Si  V  es  un  entorno  del  punto  x0,  consideraremos  en  todo  momento  f  eC(m(V) 
con  m  >  2,  una  función  convexa  en  x0.  Además,  para  obtener  una  simplificación 
en  los  cálculos,  supondremos  f(x0)  =1  ,  posteriormente  indicaremos  qué  hacer  en 
otra  situación. 

Por  otra  parte,  si  denotamos 

C((T  ir)(v)={gec(m(v)  lg(Xo)  =  r} 
y  consideramos  el  operador  normalización  a  uno  en  el  punto  x0,  dado  por 

N;c((x  o/v)  ^C(x  i)(v)  tal  que  Ntg](x)  =  9<x)  +  1  ■ 

Entonces  definimos  las  siguientes  sucesiones  de  funciones 

F-| (x)  =  l°g  f(x)  Gi(x)  =  N[F-,](x) 

F2(x)  =  logG1(x)  G2(x)  =  N[F2](x) 

Fn(x)  =  log  Gn.!  (x)  ~~  G~(x)  =  N[Fn](x) 

Así,  el  estudio  de  la  convexidad  en  el  punto  x0,  de  las  funciones  que  forman 
la  sucesión 


{Fn>n,1cC, 


(m 

(x0  ,  0) 


(V) , 


nos  va  a  permitir  medir  la  resistencia  de  la  función  f  a  ser  “concavizada” ,  median¬ 
te  sucesivas  aplicaciones  del  operador  logaritmo. 

En  la  figura  6  pueden  verse  las  gráficas  de  la  sucesión  {Fn}  para  la  función 
f(x)  =  x  -  3cos(x-1)  +  3  .  Vemos  que  F2  es  todavía  una  función  convexa  en  x0  =1 , 
Fg  tiene  en  x0  un  punto  de  inflexión  y  F4  es  cóncava  en  x0. 
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Este  índice  de  la  convexidad  va  a  establecer  lo  que  denominaremos  grado 
de  convexidad  logarítmico  de  f  en  x0  ,  grado  que  nos  va  a  indicar  cómo  es  la  con¬ 
vexidad  de  f  en  el  punto  x0  .  Para  conseguir  esto,  comenzaremos  introduciendo 
el  concepto  de  función  “n-logarítmicamente  convexa”,  concepto  que  se  obtiene 
como  generalización  del  término  de  función  logarítmicamente  convexa  estudiado 
en  el  párrafo  anterior. 


Definición  1.3. 


Diremos  que  f  es  n  -  logarítmicamente  convexa  en  x0  si  Fn  es  convexa  en  x0. 
Si  f  es  n  -  logarítmicamente  convexa  en  x0  para  todo  ne  N  diremos  que  f  es 
infinitamente  -logarítmico  convexa  en  x0. 
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A  continuación,  al  igual  que  hicimos  en  el  párrafo  anterior,  nos  interesa  carac¬ 
terizar  analíticamente  este  concepto.  Para  ello,  comenzaremos  obteniendo  una 
expresión  para  la  segunda  derivada  de  las  funciones  de  la  sucesión  {Fn}n>i  , 
derivada  que,  en  algunas  ocasiones  (carácter  estricto),  nos  permitirá  estudiar  su 
convexidad. 


Teorema  1 .4. 


Para  todo  n  e  N  se  verifica: 

(1.5)  F’n(Xo)  =  f’(xo) 

(1.6)  F”n(xo)  =  f”(x0)  -  n  [f’(x0)]2 

Demostración: 


Vamos  a  probar  ambas  expresiones  por  inducción.  Dado  que  (1.5)  nos  ser¬ 
virá  para  probar  (1 .6)  ,  demostraremos  en  primer  lugar  la  expresión  correspon¬ 
diente  a  la  primera  derivada  de  Fn.  Es  obvio  que  (1 .5)  se  verifica  para  n  =  1 ,  así 
que  supongamos  que  se  verifica  para  todo  p  e  N  con  p  <  n  ,  es  decir  F’p(x0)  = 
f’(x0).  Como  Fn(x)  =  log  Gn.-|(x)  se  sigue  que  F’n(x0)  =  G’n_i(x0)  puesto  que 


n- 1 


Jm 

^(x0  '  1 ) 


(V) 


Por  la  construcción  indicada  de  las  sucesiones  {Fn}n>i  y  {Gn}n>i  ,  es  claro  que 
G’n-l(xo)  =  F’n-l(xo)  y  aplicando  la  hipótesis  de  inducción  obtenemos  que 
F’n(x0)=G’n_i(x0)  =  f’(x0) ,  con  lo  que  (1.5)  queda  probado  . 


En  cuanto  a  (1 .6) ,  como  F-|(x)  =  log  f(x) ,  teniendo  en  cuenta  (1.1),  se  sigue 
que  F”i  (x0)  =  f”(x0)  -  [f’(x0)]2,  luego  (1 .6)  se  verifica  para  n  =  1 .  Supongamos  que 
para  todo  p  e  N  con  p  <  n  se  tiene  que  F”p(x0)  =  f”(x0)  -  p[  f’(x0)  ]2.  Entonces, 
como  Fn(x)  =  log  Gn_i(x)  =  log  [Fn_-|  (x)  +  1],  es  claro  que 


F’n(x)  = 


;,M  + 


y  F”n(x)  = 


F”n-l(X)[Fn-lW  +  1]-  F’n-l(X^ 
tFn- 1  (X>  +  1 


luego  mediante  (1 .5)  y  la  hipótesis  de  inducción  ,  como  se  verifica 


Fn--S  CC)(V)' 


obtenemos  que  F"^  =  f  "(x^  -  n  [  f  (x^  ]‘ 


lo  que  prueba  (1.6). 
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Notemos  que  si  f  es  n-logarítmicamente  convexa  en  x0  entonces  f  es  k-loga- 
rítmicamente  convexa  en  x0  para  ke  N  con  1<  k<  n. 

A  continuación,  teniendo  en  cuenta  los  diferentes  tipos  de  convexidad  que 
pueden  presentarse  en  el  punto  x0,  vamos  a  caracterizar  analíticamente  el  con¬ 
cepto  de  función  n-logarítmicamente  convexa  en  x0.  Así,  comenzaremos  anali¬ 
zando  el  caso  de  que  x0  sea  un  punto  crítico  de  f,  es  decir,  en  nuestras 
condiciones,  x0  será  un  mínimo. 


Teorema  1 .5. 

(i)  Si  xQ  es  un  mínimo  de  f  entonces  también  lo  es  de  Fn  para  todo  n  e  N. 
Además  el  carácter  del  mínimo  es  el  mismo  en  f  que  en  cualquiera  de  las  Fn,  n  e  N. 
(ii)  f  es  infinítamente-logarítmico  convexa  en  x0  sí  y  sólo  sí  x0  es  un  mínimo  de  f. 

Demostración: 

(i)  Supongamos  que  Xq  es  un  mínimo  estricto  de  f,  es  decir,  f’(xo)=0  y  f”(x0)  >  0. 
Por  (1.5)  y  (1.6)  se  sigue  que  F’n(xo)=0  y  F”n(x0)  >  0  para  cualquier  n  e  N,  por 
tanto  xQ  es  un  mínimo  estricto  de  Fn  para  todo  n  e  N. 

Si  x0  es  un  mínimo  no  estricto  de  f,  es  decir,  existe  k  e  N  número  par,  tal  que 
f’(x0)  =...  =  f(k-1(xo)  =  0yf(k(  x0)  >  0,  veamos  que  también  lo  es  de  cualquier  Fn. 
Por  (1.3),  considerando  F-|(x)  =  log  f(x),  obtenemos  que 
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Entonces  aplicando  recurrencia  a  esta  expresión  obtenemos 

n=0  f'(xo)  =  F'i(x0)  f  (xQ)  =>  F'-|(x0)  =  0 

n=1  f"(x0)  =  F'-|(x0)  f'(xo)+  F"i  (x0)  f(x0)  =>  F"-,  (x0)  =  0 

Suponi  endo  que  F^(xq)=0  para  r  <  k- 2  t  enemos  que  para 


.  0  ,  (k— 1  x  k'2fk-  2^  Jt+1  (k-2-t  .  Jk-1  ...  . 

n  =  k-2  fv  (xo)  =  t^t  V  (xo)=Fi  (xo>f<xo> 

=>  F1(k-1(xo)  =  0 
Por  otra  parte,  si  consideramos 


n  =  k  -  1 


<*o>  =  Fík(xo>,(xo>  =  Fík<xo) 


luego  F’-|(x0)  =...  =  F-^'^xo)  =  0  y  F-|(k(xo)  =  f(k(xo)  >  0.  Por  tanto  F-|  tiene  un 
mínimo  no  estricto  en  x0. 

Supongamos,  como  hipótesis  de  inducción,  que  x0  es  un  mínimo  no  estricto 
de  Fi  ,---,Fn_i .  Entonces  como  Fn(x)  =  logGn_-|(x),  por  (1.3)  se  verifica  que 


G' 


,(s+1 


n-1 


teniendo  en  cuenta  que  Gn.-|  (x)  =  Fn_i  (x)  +  1  y  que  Gn.-|  (x0)  =1  se  sigue 


(Xq)  + 


s  — 1 

I 

t  =0 


's 


(Xo)^(xo> 


y  procediendo  como  en  el  caso  n  =  1 ,  aplicando  que  x0  es  un  mínimo  no  estricto 
de  Fn-1 ,  se  obtiene  que  x0  es  un  mínimo  no  estricto  de  Fn,  para  n  e  N. 

(ii)  Si  x0  es  un  mínimo  de  f,  por  (i),  x0  es  un  mínimo  de  Fn  para  cualquier 
n  g  N,  entonces  Fn  es  convexa  en  x0  para  todo  n  e  N  y  por  tanto  f  es  infinita- 
mente-logarítmico  convexa  en  xQ. 
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Recíprocamente,  si  f  es  infinitamente-logarítmico  convexa  en  x0  entonces  Fn 
es  convexa  en  x0  para  todo  n  e  N.  De  (1.6)  se  sigue  que  f"(x0)  -  n[f'(x0)]2  >  o 
para  todo  n  e  N  y  por  (1 .5)  f'(x0)  =  F'n(x0)  =  0.  Entonces  se  sigue  que  : 

(a)  f'(x0)  =  0  y  f"(xo)>0  ó  bien 

(b)  f'(x0)  =  f"(x0)  =  0. 

El  caso  (a)  representa  la  situación  de  que  x0  sea  un  mínimo  estricto  de  f. 
Respecto  de  (b),  como  f  es  convexa  por  hipótesis,  se  obtiene  que  xQ  es  un  míni¬ 
mo  no  estricto  de  f.  # 

En  la  figura  7  puede  verse  un  ejemplo  del  caso  en  que  x0  es  un  mínimo 
estricto  de  f. 

Si  x0  no  es  un  punto  crítico  de  f,  pueden  presentarse  dos  situaciones:  f  estric¬ 
tamente  convexa  en  x0  (fig.6.)  ó  bien  f  convexa  en  sentido  no  estricto  en  x0 
(fig.8.).  A  continuación  analizamos  estas  dos  situaciones.  Denotaremos  con  E[r] 
la  parte  entera  del  número  real  r. 


Teorema  1 .6. 

Sea  f  una  función  convexa  en  sentido  estricto  en  el  punto  x0.  Entonces 

(i)  Si  r  =  f"(x0)/[  f'(x0)]2  €  N  entonces  se  verifica  que 

(a)  Para  p  =  E[r]  e  N,  tenemos  que  Fp  es  estrictamente  convexa  en  x0  y 
Fp+i  es  estrictamente  cóncava  en  x0.  Por  tanto  f  es  p-logarítmica- 
mente  convexa  en  x0  y  no  es  (p+1)-  logarítmicamente  convexa  en  x0. 

(b)  f  es  n-logarítmicamente  convexa  en  x0  para  n  <  E[r], 

(ii)  Si  r  e  N  se  verifica: 

(a)  f  es  (r-1)  -  logarítmicamente  convexa  en  x0,  siendo  Fr_i  estrictamente 
convexa  en  x0.  Además  Fr+^  es  estrictamente  cóncava  en  x0,  luego  f 
no  es  (r+1)-  logarítmicamente  convexa  en  xQ. 

(b)  f  es  r-logarítmicamente  convexa  en  x0  sí  y  sólo  si  existe  k  e  N,  núme¬ 
ro  par,  tal  que 

'=<r,_1<x0>  =  ,'<xo)t  2slsk-'  ;  F<rk_1<*o>>r<xo>k 
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En  esta  situación,  Fr  será  convexa  en  sentido  no  estricto  en  x0. 


Demostración: 

(i)  El  apartado  (a)  se  sigue  fácilmente  de  (1 .6)  ,  ya  que  si  p  =  E[r],  es  claro 
que  f"(xo)  >  p[f'(x0)]2  y  f"(x0)<(p+1)[f'(x0)]2.  En  cuanto  a  (b),  también  por  (1.6), 
resulta  obvio. 

(ii)  Al  igual  que  en  (i),  por  (1 .6)  el  punto  (a)  es  inmediato. 

Para  probar  (b),  tengamos  en  cuenta  que  Fr(x)  =  log  Gr_-|(x),  entonces  por 
(1.6)  y  como  Gr_i(x)  =  Fr_i(x)  +  1  se  sigue  que 


<1-7)  F<r^V0)=^+V0)  +  Sl(j)  F^Vo) 

t  =ov  y 


Por  otra  parte,  F"r(x0)  =  f"(x0)  -  r  [  f'(x0)]2=  0.  Entonces  f  es  r-logarítmica- 
mente  convexa  en  x0  sí  y  sólo  sí  existe  k  eN,  número  par,  tal  que  Fr”(x0)  =...  = 
F^'^Xo)  =  0  y  Fr(k(x0)  >  0,  siendo  F'r(x0)  =  f'(x0)  *  0.  De  donde  resulta  obvio 
que  Fr  es  una  función  convexa,  en  sentido  no  estricto,  en  el  punto  x0. 

Además,  por  (1 .7)  tenemos  que  : 


s=1  F"m  (x0)  =  F'r(x0)  F'r_i  (x0)  =  [  f'(x0)]2 
s=2  F"'r_1  (x0)  =  F'r(x0)  F"r.-,  (xQ)+  F"r(x0)  F'r_-,  (x0)  = 

=  F'r(x0)  F"M(xo)  =  [f'(x0)]3 
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(k-2  k-  2 

Por  recurrencia,  si  suponemos  Fr  1  (xQ)  =  [f '(xQ)]  .tenemos  que 
3  =  k  -  2  pf-’íxj  =  f  kS  (k  -  2)f  Hl(x0)Ff^2',(x0)  = 

t  =0 


k  -  2\  ík-2  k-2  k— 1 

.0  JF'r  <X0>Fr_i  (xo>  =  f'(xo)[f'(xo)]  =tf'(xo)] 


Por  otra  parte,  podemos  escribir 


s  =  k-1  ^Jx^lílx^ 

t  =0 

=  ífXo)  +  F'r  (x^-Vo)  =  F(rk(Xo)  +  [f  '(Xo)]k  >  [f '(Xo)]k  . 


Y  queda  probado  el  apartado  (b).  # 

Ejemplo  1 .2. 

f  7(0) 

El  valor  de  - r,  i  =1,2,3,  para  las  funciones 

f  ¡(0) 

x4 

f  _|(x)  =  exp  [exp(—  +  x)  -  1] 

f  2(x)  =  exp  [exp(x  -  ~)  -  1] 
x3 

f  3(x)  =  exp  [exp(—  +  x)  -  1] 

es  dos.  Si  designamos  por  {Fn¡}  la  sucesión  {Fn}  correspondiente  a  f¡,  1=1 ,2,3  se 
comprueba  (fig.9.)  que  F"2i  (0)>0,  F"22(0)  <  0  y  F23  tiene  en  el  origen  un  punto 
de  inflexión. 
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La  situación  que  nos  queda  por  analizar,  que  f  sea  una  función  convexa 
en  sentido  no  estricto  en  x0  (fig.8.),  resulta  evidente  ya  que  en  este  caso  sabe¬ 
mos  que  existe  ke  N,  número  par,  de  manera  que  f"(x0)  =....  =  f(k-1  (x0)  =  0 
y  f(k  (x0)  >  0  ;  como  además  f'(x0)  *  0,  teniendo  en  cuenta  (1 .6)  se  sigue  inme¬ 
diatamente,  que  Fn  es  una  función  estrictamente  cóncava  en  x0  para  cualquier 
neN  (fig.  8). 

Una  vez  estudiadas  las  diferentes  posibilidades  de  convexidad  que  pueden 
presentarse  en  el  punto  x0,  resulta  sencillo  observar  que  estas  quedan  perfecta¬ 
mente  determinadas  analizando  el  valor  de  la  expresión  f"(x0)[f'(x0)]"2,  que  por  (1 .6) 
nos  permite  estudiar  la  concavidad  o  convexidad  de  Fn  en  x0  para  cualquier  neN. 
Este  es  el  origen  de  la  definición  de  "grado  de  convexidad  logarítmico". 
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Definición  1.7. 

El  grado  de  convexidad  logarítmico  de  la  función  f  en  el  punto  x0  es  el  núme¬ 
ro  real  dado  por  la  expresión 


(1.8) 


cuando  x  no  es  un  mínimo 
o 

a  (m 

Notemos  que  Lf  :  C  ^  i 
convexas,  la  imagen  de  Cf 


Lf<xp)  = 


f"(xp> 

[f'(x0)]2 


de  f,  en  cuyo  caso  C  (x  )  =  +  oo. 
— >  R  ,  y  si  consideramos  funciones 

estaría  en  [0 ,  +°°  ] . 


En  el  caso  de  que  f  tome  un  valor  positivo  arbitrario  en  el  punto  x0,  conside¬ 
raremos  la  función  f*  =  f/f(xQ)  dado  que  log  f  y  log  f*  poseen  la  "misma  convexi¬ 
dad"  en  el  punto  x0,  ya  que  se  diferencian  en  una  constante  (fig.10.).  Entonces 
unos  pequeños  cálculos  nos  permiten  obtener  que  el  grado  de  convexidad  loga¬ 
rítmico  de  f  en  x0  es  el  número  real  dado  por  la  expresión 


(1.9) 


Lf(xp>  = 


<xp>  = 


f(xo)f"(xo) 

[f'(x0)]2 


si  x0  no  es  un  mínimo  de  f,  en  cuyo  caso  :  Lf(xo)  =  +°°. 


Notas: 

1 .  Si  consideramos  una  función  positiva  en  un  entorno  del  punto  x0,  es  claro 
que  el  grado  de  convexidad  logarítmico  de  una  función  convexa  es  un  número 
real  positivo,  siendo  "más  convexa"  respecto  del  logaritmo  a  medida  que  mayor 
sea  dicho  número.  Por  otra  parte,  si  bien  perdemos  el  significado  geométrico  del 
concepto,  es  claro  que  podemos  extender  formalmente  la  definición  de  Lf  al  con¬ 
junto  c(2  (V). 

2.  Es  interesante  notar  que  en  el  caso  de  que  f  no  tome  el  valor  uno  en  x0, 
la  normalización  que  hemos  considerado  no  es  la  habitual,  normalmente  se  con¬ 
sidera  N(f)(x)  =  f(x)  -  f(x0)  +1 .  Esto  se  debe  a  que  trasladar  una  función  influye  en 
la  concavidad  o  convexidad  de  su  logaritmo. 

Así,  por  ejemplo,  si  tomamos  f(x)  =  X  +  exp  (x),  x0  =  0  y  X  e  (-1 ,  +  °°),  al  con¬ 
siderar  F(x)  =  log  f(x)  obtenemos  que 
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sg(F"(x0))  =  sg  (X),  de  donde  se  sigue  que  si  trasladamos  la  gráfica  de  la  fun¬ 
ción  exponencial  de  manera  que  tome  un  valor  comprendido  entre  menos  uno  y  cero 
en  el  origen,  el  logaritmo  de  esta  función  es  una  función  cóncava  ;  siendo  convexa 
el  logaritmo  de  cualquier  traslación  de  ella  con  X  e  R+.  Por  tanto,  resulta  evidente 
que  el  grado  de  convexidad  logarítmico  no  es  una  medida  invariante  por  traslaciones. 


Como  ejemplo,  si  consideramos  f(x)  =  exp(x)  +  1  y  g(x)  =  exp(x)  -  1/2,  se 
observa  que  sus  logaritmos  poseen  distinta  convexidad  en  x0  =  0.  (fig.  1 1 .) 
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1.3.  Propiedades  del  grado  de  convexidad  logarítmico. 

En  este  apartado  vamos  a  estudiar  algunas  propiedades  del  grado  de  con¬ 
vexidad  logarítmico. 


Teorema  1 .8. 


Sea  f  eC(2(V),  positiva  y  convexa  en  V.  Entonces  se  verifica 

i)  Para  todo  número  real  r  tal  que  r  +  f(x0)  >  0, 


L 

(r 


(x„)  =  L.(xJ 
f)v  o'  fv  o' 


+  r 


xo) 


ii)  Para  todo  número  real  positivo,-  L^^(x0)  =  Lf(x0). 

iii)  Si  1/f  es  convexa  en  V  entonces  L(i /f)(x0)  =  2  -  Lf(x0). 

iv)  Sea  g  una  función  de  las  mismas  características  que  f,  de  manera  que 
esté  definida  la  composición  fog  ,  entonces 

L  (xj  =  [  Üg(xJ)  +  L  (x  )](f og)(x  ) 

(fog)  0'  1  fvav  0"  9'  0,J'  0' 

v)  Lf(x0)  >  0  y  además  Lf(x0)  =  0  sí  y  sólo  sí  existe  k  número  par,  tal  que 
f  "(x0)  =...  =  f(k_1(x0)  =  0  y  f(k(xo)>0. 

Demostración: 

Basta  calcular  las  derivadas  correspondientes  y  sustituir  en  la  expresión 
adecuada.  Para  el  apartado  (i)  tenemos 


(r+f )"(x  )(r+f )(x  ) 

L  (xj  = - - - r-5- 

<r+,)  °  [(r+f  )'(x0)]2 


f"(x0)f(x0)  f"(x0) 

[f ’(x0)]2  [f'(x0)]2 


=  Lf(xo)  +  Cf<xo)- 


¡i) 


L  (xj  = 

(rf )v  0' 


rf"(x0)rf(x0) 
[rf '(xo)]2 


f"(xo)f(xo) 

[f‘(x0)]2 


=  L,(xo)- 
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2f(xo)f'(xo)2-,"(xo),(Xo)2  1 


f(xQ) 


f(xQ) 


¡i¡)  Ll(Xo) 

f 


=  2  - 


f'(xo> 


f(x0) 


f'fc  Jf  (Xq) 

[f'íXo)]2 


=  2  -  L 


f(xo) 


iv)  L  (xj  = 
'  (fog)v  o ' 


P"(g(x0))g'(Xo)]  +  f'(g(x0))g"(x0)](fog)(Xo) 


[f  ,(g(x0))g,(x0)] 


f"(9(x0))  g"(xo) 
.|f'(g(x0))]2  +[g'(x0)]2 


(f°g)(x0)  =  [  Lf(g(x0))  +  Lg(x0)  ](fog)(x0). 


v)Lf(x0)  =  f  "(x0)  f(x0)[f'(x0)]‘2  >  0  al  ser  f  convexa  y  positiva. 

Si  Lf(x0)  =  0  se  sigue  que  f"(x0)  =  0  pues  f'(x0)  *  0  al  no  ser  x0  un  mínimo 
para  f  y  como  además  f  es  convexa  existirá  un  número  k  par  tal  que  f"(x0)  =...  = 
f(k‘1(x0)  =  0  y  f(k(x0)  >  0.  El  recíproco  resulta  obvio.  # 

Por  último  damos  unas  acotaciones  para  el  grado  de  convexidad  de  la  suma 
de  funciones. 


Teorema  1.9. 

Sean  f,  g  e  c(2(V),  positivas  y  convexas  en  V  tales  que  f'(x0)  g'(x0)  <  0,  enton- 
ces  L(f+g)(xo)  *  mín-  (  Lf(x0).  Lg(x0) }. 

Demostración: 

Podemos  suponer  Lf(x0)  <  Lg(x0)  y  estudiemos  el  signo  de  L(f+g)(x0)  -  Lf(xo). 
Es  fácil  ver  que 

(1  L  L  [gg"(f  /)2--f”f(g,)2i  +(f  f  [f  ]  -2f  f  'f  "g' 

(f+g)  f  (f’)2[f’+g’]2 
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y  como  f"(x0)  g  (Xp)  +  f(x0)  g"(x0)  >  0  por  ser  f  y  g  positivas  y  convexas  en  V,  g 
(xo)  9"(x0)  [f'(x0)r  -  f(x0)  f"(x0)  [g'(x0)]2  ^  0  por  ser  Lf(x0)  <  Lg(x0)  y  además  f 
(x0)  f'(x0)  f"(x0)  g'(x0)  <  0  por  ser  f'(x0)  g'(x0)  <  0  y  f  convexa  y  positiva,  enton¬ 
ces  L(f+g)(x0)  -  Lf(x0)  >  0  y  por  tanto  L(f+g)(x0)  >  mín.  {Lf(x0),  Lg(x0)} 


Teorema  1.10. 

Sean  f,  g  eC(2(V),  positivas  y  convexas  en  V  tales  que 
0  Lf(xo)  -  Lg(x0) 

¡i)  f(xo)  soy 

iü)  f'(xo)/  g'(x0)  <  mín.  {f(x0)/  g(x0),  f"(x0)/  g"(x0)},  entendiendo  que  si  g'(x0) 
ó  g"(x0)  son  cero,  entonces  f'(x0)  o  f"(x0)  también  lo  serán. 

Entonces  L(f+g)(x0)  <  máx.  {  Lf(x0),  Lg(x0) }  =  Lf(x0). 

Demostración: 

Basta  tener  en  cuenta  la  expresión  (1 .10)  y  que  se  verifica 

a)  g(x0)g"(xo)[f'(xo)]2  '  f(xo)f"(xo)[9'(xo)]2  *  0  ,  por  (i). 

b)  g"(x0)  f(x0)  [f'(x0)]2  -  f(x0)  f(x0)  f"(x0)g'(xo)  = 

=  f(x0)  f'(xo)  [g"(xo)  f'(xo)-f"(xo)g'(xo)]  S  0  al  ser  f(xo)f'(xo)>0 
por  (ii))  y  g"(x0)f'(xo)  -  f"(x0)g'(xo)  <  0  Por  ("')• 

c)  f"(xo)  g(x0)  [f'(xo)]2  '  W  f'(xo)  f"(xo)9'(xo)  - 

=  f(x0)  f"(x0)  [  g(xQ)  f(x0)  -  f(x0)g'(x0)]  <  0  por  razones  análogas  a  las 
utilizadas  anteriormente. 

De  a),  b)  y  c)  se  deduce  que  L(f+g)(x0)  -  Lf(x0)  <0.  # 
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CAPÍTULO  2 

MEDIDA  DE  LA  CONCAVIDAD  DE  UNA  CURVA 
MEDIANTE  EL  OPERADOR  EXPONENCIAL 

2.0.  Introducción 

Se  puede  realizar  un  análisis  análogo  al  del  Capítulo  1  para  estudiar  la  con¬ 
cavidad  de  una  función,  utilizando  la  acción  del  operador  exponencial  sobre  las 
funciones  cóncavas.  Al  ser  la  exponencial  una  función  convexa,  va  disminuyendo 
la  concavidad  de  la  función  a  la  que  se  le  aplica  hasta  conseguir,  que  en  una  de 
sus  repetidas  aplicaciones,  la  función  obtenida  sea  convexa.  Obtendremos  así  un 
índice  de  medida,  en  cada  punto,  de  la  concavidad  de  la  función  considerada  . 


La  figura  12  muestra  que  hemos  de  aplicar  tres  veces  el  operador  exponencial 
a  la  función  f(x)  =  x  -  x^,  que  es  cóncava  en  el  punto  x0  =0,  para  que  deje  de  serlo  . 
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2.1.  Funciones  exponencialmente  cóncavas 

Sea  V  un  entorno  del  punto  x0  y  f  e  c(m(V)  con  m  >  2,  cóncava  en  V.  Deci¬ 
mos  que  f  es  una  función  exponencialmente  cóncava  en  x0  si  exp  f  es  una  fun¬ 
ción  cóncava  en  x0.  Entonces,  si  f  es  una  función  exponencialmente  cóncava  en 
x0,  existe  un  entorno  W  de  x0  y  una  función  h  e  C¡(m(W)  cóncava,  W  c  V,  de 
manera  que  f(x)  =  log  h(x)  en  W,  es  decir,  h(x)=e*(x)  en  W.  Operando  con  la 
expresión  obtenida  para  la  función  f  se  sigue  que  h'(x)  =  f'(x)  exp  [f(x)]  =  h(x)  f'(x) 
y. además  h"(x)  =  f"(x)  h(x)  +  [f'(x)]2  h(x)  =  [f"(x)  +  f'(x)2]  h(x).  De  estas  igualda¬ 
des  se  obtiene 

(2.1)  h"(x)  =  [f"(x)  +  f(x)2]  ef(x) 

igualdad  que  nos  va  a  permitir  realizar  el  estudio  de  la  concavidad  de  la  función  h. 
Obtenemos  así  los  siguientes  resultados,  análogos  a  los  Teoremas  1.1  y  1.2. 


Teorema  2.1. 

(i)  Si  x0  no  es  un  punto  crítico  de  f,  se  verifica 

(a)  f  es  exponencialmente  cóncava  en  x0,  en  sentido  estricto, 


sí  y  sólo  sí 


nv 

f,(x0)2 


>  1. 


(b)  f  es  exponencialmente  cóncava  en  x0,  en  sentido  no  estricto  sí  y  sólo 
sí  existe  keN  número  par,  k  <  m,  de  manera  que 


(2.2) 


|f(t(x0)=(-1)t+1(t-1)![f'(x0)]t 
f(k(x0)<-(k-1)l  [f'(Xo)]k 


2  <  t  <  k  -  1 


(ii)  Si  f  es  una  función  cóncava  en  sentido  no  estricto,  entonces  f  no  es  expo¬ 
nencialmente  cóncava. 


Demostración: 


Como  h'(x0)  =  f'(x0)  exp  [f(x0)]  *  0,  xQ  no  es  punto  crítico  de  h.  El  apartado 
(a)  del  punto  (i)  es  claro,  por  (2.1).  Para  probar  (b),  como  h  =  ef(x)  es  cóncava  en 
x0,  pero  en  sentido  no  estricto,  se  verifica  que  existe  keN  número  par,  k  <  m,  de 
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manera  que  h"(x0)  =...  =  h^‘"*(x0)  =  0  y  h(^(x0)  >  0,  luego  teniendo  en  cuenta 
(2.1)  se  obtiene  que  -  f"(x0)=[  f'(x0)]^. 

Para  probar  (2.2),  basta  aplicar  un  proceso  inductivo,  teniendo  en  cuenta  que 
h'(x)  =  h(x)  f'(x) ,  por  tanto  se  verifica 

(2.3)  hv  (xQ)=  I  L  M  (X0)I'|V  (Xq). 

t  =0  ■  ' 

El  recíproco  se  obtiene  fácilmente  a  partir  de  (2.3). 

La  primera  de  las  condiciones  (2.2)  no  es  suficiente  para  que  la  función  f  sea 
exponencialmente  cóncava,  como  prueba  el  ejemplo  siguiente. 

Ejemplo  2.1. 

Las  funciones  f  (x)  =  Log .  (1  -  x  +  *-) 
f2(x)  =  Log.(1-x-^-) 
f  3(x)  =  Log.(1-  x-  ^-) 

verifican  dicha  condición  para  x0  =  0  y  k  =  2.  Sin  embargo  exp  f-j  es  convexa,  exp 
f2  es  cóncava  y  exp  tiene  en  x0  un  punto  de  inflexión. 
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Teorema  2.2. 


Si  xQ  es  un  punto  crítico  de  f,  es  decir  un  máximo,  entonces  se  verifica 


x_>x°  [f  ’(x)] 

Además  xQ  será  un  máximo  de  h  =  e*(x),  manteniéndose  para  h  el  carácter 
estricto  o  no  estricto  que  x0  tenga  para  f. 

Demostración: 


Para  probar  este  teorema  es  necesario  distinguir  dos  situaciones  posibles.  A 
parte  de  verificarse  que  f'(x0)  =  0  puede  ocurrir  que 

1  ■)  f"(x0)  <  0 

2.)  Exista  keN  número  par,  k  <  m,  de  manera  que 
f(x0)=...  =  f(k-i(x0)  =  oy  f(k(x0)<o. 

En  el  caso  1.)  resulta  obvio  que,  como  f"(x0)  <  0  se  verifica  (2.4).  Además, 
se  tiene  que  h'(x0)  =  0  y  h"(x0)  <  0,  luego  se  mantiene  el  carácter  estricto  del 
máximo  x0  de  f  para  la  función  h  =ef(x). 

En  cuanto  a  2.),  aplicando  la  regla  de  L'Hópital,  se  prueba  que  también  veri¬ 
fica  (2.4).  Y  considerando  la  expresión  (2.3)  se  obtiene  que  h'(x0)  =  h"(x0)  =...  = 
h(k_1  (xD)  =  0  y  h(k(x0)  <  0,  luego  x0  es  un  máximo  no  estricto  para  h.  # 

Notemos  que  del  punto  (a),  del  apartado  (i)  de  este  resultado  se  sigue  que  f 
ha  de  ser  estrictamente  cóncava  en  x0.  Por  otra  parte,  es  sencillo  probar  que  si  f 
es  cóncava  en  xQ  de  forma  no  estricta,  no  siendo  x0  un  punto  crítico  de  f,  enton¬ 
ces  h  =  e^x)  es  una  función  convexa  en  x0.  Para  demostrar  esto  basta  conside¬ 
rar  la  expresión  (2.1). 
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2.2.  Definición  del  grado  de  concavidad  exponencial  de  una 
función 


Para  dar  un  índice  de  medida  puntual  de  la  concavidad  de  una  función,  repe¬ 
timos  el  proceso  realizado  en  el  caso  de  la  convexidad  utilizando  el  operador 
exponencial  en  vez  del  logarítmico. 

Sea  V  un  entorno  del  punto  x0  y  f  e  c(  m(V)  con  m  >  2,  una  función  cónca¬ 
va  en  x0.  Además,  para  obtener  una  simplificación  en  los  cálculos,  suponemos 
f(xQ)  =  0,  posteriormente  indicaremos  qué  hacer  para  un  valor  arbitrario  de  f  en  el 
punto  x0. 

Llamamos  T  al  operador  normalización  a  cero,  dado  por 


T  :C(x  i)^  "*C(X  0)^  tal  qUS  T[gl(x)  =9M  -9(xo) 


Entonces  definimos  las  siguientes  sucesiones  de  funciones: 


Q-l  (x)  =  exp  f(x) 
Q2(x)  =  exp  Ri  (x) 


Rl(x)=T[Q1](x) 
R2(x)  =  T[Q2](x) 


Qn(x)  =  exp  Rn .  i  (x) 


Rn(x)  =  T[Qn](x) 


Así,  el  estudio  de  la  concavidad  en  el  punto  x0,  de  las  funciones  que  forman 
la  sucesión 

|Qn>nacC0,1)OT- 


nos  va  a  permitir  medir  la  resistencia  de  la  función  f  a  transformarse  en  convexa, 
mediante  sucesivas  aplicaciones  del  operador  exponencial.  Así  por  ejemplo,  si 
consideramos  la  función  f(x)  =  sen(x-l)  -  x^+i  en  el  punto  x0  =1,  se  comprueba 
que  Qi  y  Q2  son  cóncavas  en  x0  mientras  que  Q2  ya  pasa  a  ser  convexa  en 
el  punto  x0.  (fig.  14.) 
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Notemos  que 

'  o  ’  ' 

Esta  medida  nos  permite  definir  el  grado  de  concavidad  exponencial  de  f  en 
x0.  Introducimos  para  ello  el  concepto  de  función  “n-exponencialmente  cóncava” 


Definición  2.2. 

Diremos  que  f  es  n  -  exponencialmente  cóncava  en  x0  si  Qn  es  cóncava  en  x0. 

Si  f  es  n  -  exponencialmente  cóncava  en  x0  para  todo  neN  diremos  que  f 
es  infinitamente  -  exponencialmente  cóncava  en  x0. 
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Tratamos  ahora  de  caracterizar  analíticamente  este  concepto.  Como  en  el 
caso  de  la  convexidad  comenzaremos  obteniendo  una  expresión  para  la  segun¬ 
da  derivada  de  las  funciones  de  la  sucesión  {Qn}n>-| . 


Teorema  2.3. 

Para  todo  ne  N  se  verifica: 

(2.5)  Q'n(x0)  =  f’(x0) 

(2.6)  Q”n  (Xq)  =  f”(x0)  +  n  [f’(x0)]2 
Demostración: 

Se  demuestra  por  inducción,  de  manera  análoga  al  Teoremal  .4.  # 

Notar  que  si  f  es  n-exponencialmente  cóncava  en  x0  entonces  f  es  k-expo- 
nencialmente  cóncava  en  x0  para  keN  con  1<  k<  n. 


A  continuación  caracterizamos  analíticamente  el  concepto  de  función  n-expo¬ 
nencialmente  cóncava  en  x0. Distinguimos  dos  situaciones,  en  primer  lugar  anali- 
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zaremos  el  caso  en  que  x0  sea  un  punto  crítico  de  f,  es  decir,  un  máximo  (fig.15) 
y  posteriormente  el  caso  en  que  x0  no  sea  un  punto  crítico  ( fig.  14  y  16  ). 


Teorema  2.4. 

(i)  Si  x0  es  un  máximo  de  f  entonces  también  lo  es  de  Qn  para  todo  ne  N.  Ade¬ 
más  el  carácter  del  máximo  es  el  mismo  en  f  que  en  cualquiera  de  las  Qn,  ne  N. 

(ii)  f  es  infinítamente-exponencialmente  cóncava  en  x0  sí  y  sólo  sí  x0  es  un 
máximo  de  f. 

Demostración: 

Es  análoga  a  la  del  Teorema  1 .5.  # 


Como  en  el  caso  de  la  convexidad,  cuando  x0  no  es  un  punto  crítico  de  f  se 
nos  presentan  dos  posibles  situaciones  :  f  estrictamente  cóncava  en  x0  ( fig.  14.) 
ó  bien  f  cóncava  en  sentido  no  estricto  en  x0  (fig. 16.). 

El  Teorema  siguiente  recoge  los  resultados  obtenidos  en  estas  situaciones, 
cuya  demostración  es  análoga  a  la  del  Teorema  1 .6. 
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Teorema  2.5. 

Supongamos  que  x0  no  es  un  punto  crítico  de  f,  que  es  una  función  cóncava 
en  sentido  estricto  en  el  punto  x0.  Entonces 

(i)  Si  r  =  -  f”(x0)/[  f’(x0)]2  g  N  entonces  se  verifica  que 

(a)  Sea  p  =  E[r]  e  N,  entonces  Qp  es  estrictamente  cóncava  en  xQ  y  Qp+1 
es  estrictamente  convexa  en  xQ  Por  tanto  f  es  p  -  exponencialmente 
cóncava  en  xQ  y  no  es  (p+1)  -exponencialmente  cóncava  en  xQ. 

(b)  f  es  n-exponencialmente  cóncava  en  x0  para  n  <  E[r]. 

(ii)  Si  r  e  N  se  verifica: 

(a)  f  es  (r-l)-exponencialmente  cóncava  en  xQ,  siendo  Qr_i  estrictamente 
cóncava  en  x0.  Además  Qr+i  es  estrictamente  convexa  en  x0,  luego 
f  no  es  (r+1)  -  exponencialmente  cóncava  en  x0. 

(b)  f  es  r  -  exponencialmente  cóncava  en  x0  si  y  sólo  si  existe  k  €  N  , 
número  par,  tal  que 

(t  (t  t+1  t 

Qr-i(xo>=f  (x0)  =  (-1)  (t  -  l)lf'(Xo)  2  <  t  <  k  -  1 

Q(rk_l(xo)<-(k“1)!f'(xo)k-  # 

En  la  figura  13  puede  verse  la  situación  del  apartado  (ii). 

La  situación  que  nos  queda  por  analizar,  que  f  sea  una  función  cóncava  en 
sentido  no  estricto  en  x0,  resulta  evidente  ya  que  en  este  caso  sabemos  que  exis¬ 
te  keN,  número  par,  de  manera  que  f”(x0)  =....  =  f ( k " 1  (x0)  =  0  y  f(k  (x0)  <  0  ; 
como  además  f’(xo)*0,  teniendo  en  cuenta  (2.6)  se  sigue  inmediatamente  que  Qn 
es  una  función  estrictamente  convexa  en  x0  para  cualquier  neN.  (fig.  16.) 

Una  vez  estudiadas  las  diferentes  posibilidades  de  concavidad  que  pueden 
presentarse  en  el  punto  x0,  resulta  sencillo  observar  que  estas  quedan  perfec¬ 
tamente  determinadas  analizando  el  valor  de  la  expresión  -  f”(x0)[f’(x0)]'  2.  la 
cual  por  (2.6),  nos  permite  estudiar  la  concavidad  o  convexidad  de  Qn  en  xQ 
para  cualquier  neN.  Este  es  el  origen  de  la  definición  de  “grado  de  concavidad 
exponencial”. 
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Definición  2.6. 

El  grado  de  concavidad  exponencial  de  la  función  f  en  el  punto  x0  es  el  núme¬ 
ro  real  dado  por  la  expresión 


(2.7) 


f"(xo) 

[f'(x0)]2 


cuando  x0  no  es  un  máximo  de  f,  en  cuyo  caso  Éf(xQ)  =  +°o . 


A  _(m 

Notemos  que  E  :  ^  R  ,  si  consideramos  funciones  cóncavas, 

la  imagen  estaría  en  [0,  +°°]- 

En  el  caso  de  que  f  tome  un  valor  arbitrario  en  el  punto  x0,  consideraremos 
la  función  f*(x)  =  f(x)  -  f(x0).  Dado  que  exp  f  y  exp  f*  poseen  la  “misma  concavi¬ 
dad”  en  el  punto  xQ  ,  ya  que  exp  f*(x)  =  [exp  f(xQ)  ]"  "'exp  f(x),  y  por  tanto  se  dife¬ 
rencian  en  el  producto  por  una  constante  positiva.  Entonces  el  grado  de 
concavidad  exponencial  de  f  en  x0  es  el  número  real  dado  por  la  expresión 

A  f"(X0) 

(2.8)  E,(x0)  =  É,,  (x0)  = - ¡4 

|f(x0)] 


si  x0  no  es  un  máximo  de  f,  en  cuyo  caso  Ef(x0)  =  +°°. 

Notemos  que  si  consideramos  una  función  en  un  entorno  del  punto  xQ,  es 
claro  que  el  grado  de  concavidad  exponencial  de  una  función  cóncava  es  un 
número  real  positivo.  Siendo  “más  cóncava”  respecto  a  la  función  exponencial  a 
medida  que  mayor  sea  dicho  número. 

Por  otra  parte,  resulta  evidente  que  el  grado  de  concavidad  exponencial  sí 
es  invariante  por  traslaciones,  propiedad  geométrica  que  nos  va  a  ser  de  gran 
utilidad. 
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2.3.  Propiedades  del  grado  de  concavidad  exponencial. 

Análogamente  a  como  se  probaron  las  propiedades  del  grado  de  convexidad 
logarítmico  en  el  párrafo  1.3  del  Tema  I,  unos  sencillos  cálculos  nos  dan  las 
siguientes  propiedades  para  el  grado  de  concavidad  exponencial. 


Teorema  2.7. 

Sea  f  ed  2(V)  cóncava  en  V,  entonces  se  verifican  las  siguientes  propiedades 

i)  Para  todo  número  real  r,  E(  r  +  f  )(x0)  =  Ef  (x0). 

ii)  Para  todo  número  real  positivo,  E(  rf  )(x0)  =  Ef  (x0)/r. 

iii)  Si  1/f  es  cóncava  en  V  ,  E(  i  /f  )(x0)  =  f(xQ)[  Lf  (x0)  -  2  ]  = 

=  -  f(x0)[  2  +  f(x0)Ef  (x0)] 

iv)  Ef(x0)  >  0  y  además  Ef(x0)  =  0  sí  y  sólo  sí  existe  un  número  par  k  tal  que 
f”(x0)  =...  =  f( k " 1  (x0)  =  0  y  f(  k(x0)  <  0. 

v)  Sea  g  una  función  de  las  mismas  características  que  f  de  manera  que  esté 
definida  la  composición  fog  ,  entonces 

E(  fog)(xo)  =  Ef  (9(xo))  +  Eg(x0).  # 

Como  ocurría  con  el  grado  de  convexidad,  se  tienen  las  siguientes  acotacio¬ 
nes  para  el  grado  de  concavidad  de  la  suma 


Teorema  2.8. 

Sean  f,  g  eC(2(V),  cóncavas  en  V  tales  que  f’(xo)g’(xo)<0,  entonces 
E(f+g)(xo)  *  mín-  íEf(xo).  Eg(x0)  }■ 

Demostración: 

Se  supone  Ef  (x0)  <  Eg(x0)  y  se  estudia  el  signo  de  E(f +g)(x0)-Ef(x0).  # 

Teorema  2.9. 

Sean  f,  geC(2(V),  cóncavas  en  V  tales  que 

0  Ef  (xo)  -  Eg(xo) 

¡i)  E(xo)  9'(xo) 

Entonces  E(f+g)(x0)  <  máx.  {Ef  (x0),  Eg(x0) }  =  Ef  (x0).  # 
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CAPÍTULO  3 

ÍNDICES  DE  MEDIDA  DE  LA  CONVEXIDAD 
Y  DE  LA  CONCAVIDAD 


3.0.  Introducción. 

En  este  capítulo  vamos  a  considerar  una  función  cóncava  (p  en  un  dominio 
W  c  R,  <p  e  c(P(W),  y  f  e  c(P(V) ,  p  >  2,  una  función  convexa  en  un  entorno  V 
de  un  punto  x0,  con  y0  =  f(  x0)  e  W. 

La  curvatura  de  una  función  f  en  x0  [14],  que  viene  dada  por 

f"(Xo) 

la  expresión  K(f)(xJ  = - „  es  una  medida  de  la 

[l+f'ÍXo)  ] 

convexidad  de  una  función  convexa  f  en  cada  punto.  Es  claro  que  si  ip'(yo)=1 
entonces  cpof  tiene  curvatura  más  pequeña  que  f,  puesto  que 

(3.1.)  (cpof)'(x0)  =  tp'(y0)  f  '(x0)  y 

(3.2.)  ( <pof)"(x0)  =  <p"(y0)  f'(x0)2  +  cp'(y0)f"(x0). 

( <pof  )  "(X  ) 

As.  K((pof)(xo)  =  - —  <  K(f)(xo)  . 

[1  +(<P°f)'(x0)  ] 

Por  tanto  aplicando  un  operador  cóncavo  a  una  función  convexa  obtenemos 
una  función  de  curvatura  más  pequeña  que  la  función  convexa  inicial.  Podemos 
observar  esta  situación  en  el  siguiente  ejemplo 
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Ejemplo  3.1. 

Sea  f(x)  =  ex  -1  y  cp(x)  =  1  +  x  -  x2/4,  entonces  cp  o  f(x)  =  ex  -  (ex  - 1)2/4.  La 
fig.  17.  compara  K(f)  y  K((pof)  para  x0  =  0. 

En  este  estudio  consideramos  el  efecto  del  operador  cóncavo  <p  sobre  el  con¬ 
junto  de  las  funciones  convexas  reales  en  un  entorno  de  un  punto  dado.  Es  decir, 
de  qué  manera  el  operador  <p  “concaviza”  a  la  función  f  a  la  que  se  le  aplica. 


Introducimos  un  índice  de  medida  de  la  convexidad  para  cada  función  con¬ 
vexa  f  de  c(P(V)  en  un  punto  dado.  Este  índice  denotará  el  número  de  veces  que 
hemos  de  componer  la  función  f  con  el  operador  <p  para  obtener  una  función 
cóncava. 

Terminamos  el  capítulo  estudiando  de  forma  análoga  la  acción  de  un  opera¬ 
dor  convexo  sobre  el  conjunto  de  las  funciones  cóncavas  de  c(P(V). 
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3.1.  Funciones  cp  -  Convexas. 

Con  la  notación  anterior  diremos  que  f  es  (p-convexa  en  xQ  sí  y  sólo  sí  tpof 
es  una  función  convexa  en  x0.  Como  la  convexidad  de  ipof  puede  ser  en  sentido 
estricto  o  no  estricto,  es  decir,  que  (ípof)”(Xg)  >  0  o  bien  que  exista  un  número 
par  k  e  N  tal  que  (cpof)”(x0)  =...  =  (cp0f)C  ^ _  '(x0)  =  0  y  (ipof)(  ^(x0)  >  0,  tenemos 
así  los  conceptos  de  función  (p-convexa  en  sentido  estricto  y  no  estricto.  En  lo  que 
sigue  también  utilizaremos  otros  conceptos  como  el  de  punto  mínimo  x0  en  sen¬ 
tido  estricto  y  no  estricto.  Su  definición  es  análoga,  es  decir,  si  f’(x0)  =  0  y 
f”(x0)  >  0  o  bien  f’(x0)  =  f”(x0)  =...  =  f(  (x0)  =  0  y  f(  ^(x0)  >  0  donde  k  es  un 
número  par. 

En  primer  lugar  queremos  caracterizar  las  funciones  <p-convexas.  Notar  que 
si  (p  es  una  función  decreciente  en  un  entorno  de  y0,  cp‘(yo)  -  0.  entonces  se  sigue 
de  (3.2.)  que  (pof  será  siempre  una  función  cóncava.  Se  obtiene  el  mismo  resul¬ 
tado  si  cp  es  una  función  cóncava  en  sentido  no  estricto  en  y0  =  f(x0),  pues  enton¬ 
ces  <p"(y0)  =0.  Para  no  perder  el  sentido  geométrico  inicial,  en  lo  que  sigue 
supondremos  que  cp  es  una  función  creciente  y  cóncava  en  sentido  estricto  en  un 
entorno  W  de  y0,  es  decir,  (p‘(y)>0  y  <p"(y)  <  0  en  un  entorno  W  de  y0. 

Considerando  el  valor  de  las  expresiones  (3.1.)  y  (3.2.)  en  x0  y  los  diferen¬ 
tes  tipos  de  convexidad  obtenemos  los  siguientes  resultados 


Teorema  3.1. 


Sea  xQ  tal  que  f’(x0)  *  0. 

(i)  Si  f  es  una  función  estrictamente  convexa  en  x0  se  verifica 

(a)  f  es  strictamente  cp-  convexa  en  xQ  sí  y  sólo  sí 

<p'(y0)  f"(x  ) 

- 5 - Q— r-  >1. 

[-  <P"(y0)]f'(x0) 

(b)  Si  f  es  una  función  cp  -  convexa  en  x0  en  sentido  no  estricto,  entonces 


(3.3) 


<p'(y0)  f"(x0) 

[-  <p"(y0)]  f'(x0)2 
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(i¡)  Si  f  es  una  función  convexa  en  x0  en  sentido  no  estricto,  entonces  f  no  es 
9  -  convexa. 

Demostración: 

(a)  y  (b)  se  siguen  fácilmente  de  (3.1.)  y  (3.2.). 

Al  ser  f”(xo)=0  implica  que  ((pof)”(x0)  =  cp”(y0)  f’(x0)2  <  0  y  se  sigue  ii).  # 


Nota: 

La  condición  (3.3.)  no  es  suficiente  para  que  la  función  f  sea  cp-convexa  en 
sentido  no  estricto,  como  se  puede  ver  en  el  siguiente  ejemplo. 


Ejemplo  3.2. 

Dadas  las  funciones  f¡  ,  i  =  1 ,2,3,  por  las  expresiones 

f1(x)=exp(^-  +  x),  f2(x)  =exp(x-^-)  y  f  3(x)  =  exp  ( +  x)  . 

Si  tomamos  como  (p(x)  =  log  x  ,  estas  funciones  verifican  la  condición  (3.3.) 
para  xo=0.  Denotemos  F¡  =  cp  of¡ ,  i  =  1 , 2,  3. 

La  f i g .1 8 .  muestra  que  F^(0)  >  0  ,  F"(0)  <0  y  F^  tiene 
un  punto  de  inflexión  en  x0. 
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Por  otra  parte,  si  f  tiene  un  mínimo  en  x0  entonces  tenemos 

q>'(f(x))f"(x) 


Ko  [-  cp"(f(x))]f'(x) 


Además,  cpof  tendrá  un  mínimo  en  x0  manteniéndose  el  carácter  estricto  o  no 
estricto  que  f  tiene  en  x0. 


53 


LA  CONVEXIDAD  EN  LA  RESOLUCIÓN  DE  ECUACIONES  ESCALARES  NO  LINEALES 


3.2.  (p-grado  de  convexidad. 

En  este  párrafo  definiremos  un  índice  puntual  de  medida  de  la  convexidad  de 
una  función.  Como  ya  comentamos  en  la  introducción  aplicaremos  reiteradamen¬ 
te  el  operador  9  a  la  función  f  y  definiremos  este  índice  como  el  número  de  veces 
que  necesitamos  aplicar  dicho  operador  hasta  obtener  una  función  cóncava. 

En  las  condiciones  usuales,  definimos  las  siguientes  sucesiones: 

F 1  (x)  =  cp  of(x)  G-)  (x)  =  N[  F-j  ](x)  =  9  of(x)  -  9  of(x0)  +  f(x0) 

F2(x)  =  9  oGt  (x)  G2(x)  =  N[  F2](x)  =  F2(x)  -  F2(x0)  +  f(x0) 


Fn(x)  =  9oGn.1(x)  Gn(x)  =  N[  F n](x). 

Notemos  que  Fk(x0)  =  9of(x0)  y  Gk(x0)  =  f(x0)  para  cada  keN.  Así,  estudia¬ 
remos  la  convexidad  de  las  funciones  Fn  para  medir  la  resistencia  que  tiene  la 
función  f  a  ser  “concavizada”  por  el  operador  9.  Para  realizar  este  estudio  intro¬ 
duciremos  el  concepto  de  función  n-9-convexa,  que  es  una  generalización  del 
concepto  de  función  9-convexa  que  hemos  estudiado  en  el  último  párrafo. 


Definición  3.2. 

f  es  una  función  n  -  9  -  convexa  en  x0  si  Fn  es  convexa  en  x0. 

Si  f  es  una  función  n  -  9  -  convexa  en  x0  para  todo  ne  N  ,  diremos  que  f  es 
indefinidamente  9-convexa. 

A  continuación  vamos  a  caracterizar  analíticamente  estos  conceptos  obte¬ 
niendo  una  expresión  de  las  derivadas  sucesivas  de  Fn,  para  todo  ne  N. 


Teorema  3.3. 

Se  verifica  para  todo  n  e  N 
0)  F'n  (x0)  =  9'(y0)n  f'(x0). 

(i¡)  Fn(xo)  =  (p^ 1  (P'(y0)V'(x0)2+  9'(y0)nf"(x0) 

k=0 
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Demostración: 

Aplicaremos  inducción  para  probar  ambas  igualdades.  Las  dos  fórmulas  son 
claras  para  n  =1 .  Supongamos,  por  hipótesis  de  inducción,  que  (i)  y  (ii)  se  verifi¬ 
can  para  m.  Como  F’m  +  i(x)  =  (p‘(Gm(x))  F’m(x),  aplicando  la  hipótesis  de  induc¬ 
ción  se  sigue  que  F’m  +  i(x0)  =  cp‘(y0)m  + V(x0).  Esto  prueba  (i). 

Por  otra  parte,  al  ser  F”m  +  1(x)=(p”(Gm(x))Fm(x)2+cp’(Gm(x))F”m(x),  un 
argumento  de  inducción  prueba  (ii).  # 

Notemos  que  si  f  es  n-cp-convexa  en  x0,  entonces  f  es  k-cp-convexa  en  x0  para 
todo  1  <  k  <  n. 

Cuando  Fn  sea  convexa  en  sentido  no  estricto  calcularemos  las  derivadas 
sucesivas  por  recurrencia  aplicando  (3.1 .)  y  (3.2.),  en  otro  caso  tendremos 


(3.4.) 


F(nk+V)=  I 

t=0 


[^(G^x))] 


(t 


F(k-1 +1 

n-1 


A  continuación,  teniendo  en  cuenta  los  diferentes  tipos  de  convexidad  en  x0, 
caracterizaremos  el  concepto  de  función  n-(p-convexa  en  x0  Comenzaremos 
estudiando  el  caso  en  que  f  tenga  un  mínimo  en  xQ. 


Teorema  3.4. 


(i)  Si  f  tiene  un  mínimo  en  x0,  Fn  también  lo  tiene  para  todo  neN.  Además  el 
carácter  del  mínimo  es  el  mismo  para  f  y  para  cada  Fn,  neN. 

(ii)  f  tiene  un  mínimo  en  x0  sí  y  sólo  sí  f  es  indefinidamente  cp-convexa  en  x0. 
Demostración: 


(i)  es  inmediato  si  x0,  es  un  mínimo  estricto.  En  otro  caso,  es  decir,  si  existe 
un  número  par  k  e  N  tal  que  f’(x0)  =..  =  f(  k-  ^  (x0)  =0  y  f(  ^(x0)  >  0,  aplicaremos 
inducción  para  probar  que  Fn  tiene  un  mínimo  en  x0  para  todo  n  e  N,  en  sentido 
no  estricto.  Para  n  =1 ,  como  F’i  (x)  =  tp‘(f(x))f’(x),  se  tiene  que 


(r  +  1  ir]  (t  (r-t  +  1  (k-1 

;  (x)  =  2,  t  wwx))]  f  (x)  y  f’1(x0)  =  fiii(x0)  =  ...  =  f;  (x0)= 

t  =0  ^  ' 


=  0.  por  tanto  =  cp'(y0)  f  (k(xQ)  >  0. 
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Obtenemos  así  que  F-|  tiene  un  mínimo  en  x0  en  sentido  no  estricto.  Un  pro¬ 
ceso  inductivo  nos  asegura  que  Fn  tiene  un  mínimo  en  x0  en  sentido  no  estricto 
para  todo  n  e  N. 

(ii)  se  sigue  fácilmente  a  partir  de  (i).  # 

Cuando  x0  no  es  un  punto  crítico  de  f,  pueden  ocurrir  dos  casos:  f  es  estric¬ 
tamente  convexa  en  x0  o  f  es  convexa  en  xQ  en  sentido  no  estricto. 


Teorema  3.5. 


Sea  f  una  función  estrictamente  convexa  en  x0  con  f’(xo)*0,  así  se  tiene 
(i)  Si  se  verifica  que 

> 

(t>: 

entonces  Fn  es  estrictamente  convexa  en  x0  y  Fn+-|  es  estrictamente 
cóncava  en  xQ.  Por  tanto  f  es  n  -  <p  -  convexa  y  no  es  (n+1 )  -  cp  -  convexa. 


n 

5>(y0: 


k=0 


)  >  [- 
(a) 


<P'(yo)  f  "(x0) 

<p"(y0) f  '(xo)£ 


w 


<p'(y0)  f"(x0)  ^  k 

(Í  Í)  Si  [  -  ] - r=Z^(P'(yo)  entonces  sf  es  (n- 1 )-  <p- 

*1*0*  f'(x  )  ° 


-  convexa  en  x 


o- 


Además  Fn .  -j  es  estrictamente  convexa  en  xQ,  y  Fn+ 1  es  estrictamente  cón¬ 


cava  en  xr 


Demostración: 

De  la  parte  (ii)  del  Teorema  3.3.,  se  sigue  fácilmente  que 


Fn(*0> 


n- 1 


<p"(yJ<p'(yJn  V(xn)2  k=o 


=  S(p,(y0)' 


<p'(y0)f'(xo) 

9"(y0)f'(x0)2 


<  0,  por  (b) 
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y  como  el  denominador  es  menor  que  cero  se  sigue  que  F”n(x0)  >  0  y  por  tanto 
Fn  es  una  función  estrictamente  convexa. 

Análogamente  al  caso  anterior,  se  comprueba  fácilmente  que 


■F«1<X0> 


<p"(yn)<p'(y  /nxj2 


X  <p,<y0)k  + 

k=0 


<P'(y0)f"(x0) 

<p"(y0)  f'(x0)2 


>o  , 


por 


(a) 


y  como  F”n+1(x0)  <  0,  entonces  Fn+i  es  una  función  estrictamente  cóncava. 
Esto  prueba  (i). 

Para  probar  (ii)  deducimos  de  las  hipótesis  que 


2n— ^  ? 

Fn-i(xo)=-(p"(y0)(p'(y0)  n*o>  >0 


y  por  ser  cp  "(y q)  <  0  , 


obtenemos  que  Fn .  ■)  es  estrictamente  convexa. 

C°mo  Fn(xo)  =  °  y  F”n+l(xo)  =  <P"(yo>  V(y<fn*  <0  se 

sigue  que  Fn  +  i  es  una  función  estrictamente  cóncava.  # 

Este  último  resultado  nos  permite  afirmar  que  la  cantidad 

<p'(y0)  f"(x0)  ^  k 

l - 1  y  su  relacion  con  s  =  2ii  ^'(Vn)  ’  establece  el 

tp"(y0)  f'(x0)  k=o 

9  -  grado  de  convexidad  de  f  en  x0.  Como  este  viene  definido  por  el  número  natu¬ 
ral  n  tal  que  Fn  sea  convexa  y  Fn+i  no  lo  sea,  teniendo  en  cuenta  los  diferentes 
valores  de  s  parece  natural  normalizar  el  operador  cp  de  manera  que  cp'(y0)  =  1 . 
Esta  normalización  no  cambia  las  propiedades  geométricas  de  cp,  es  decir,  cp  es 
creciente  y  cóncava  en  un  entorno  W  de  y0. 

Por  otra  parte,  en  esta  situación,  nuestra  interpretación  geométrica  del 
cp-grado  de  convexidad  en  relación  con  la  curvatura  es  clara.  Dado  el  operador  cp 
en  las  condiciones  indicadas,  consideremos  vp(y)  =  cp(y)/cp’(y0)  y  repetimos  el  pro¬ 
ceso  anterior  para  este  operador  \p,  entonces  calculando  las  correspondientes 
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sucesiones  Hn(x)  =  yoKn .  i  (x)  donde  Kn(x)  =  N[Hn](x),  H-j  (x)  =  y  o  f(x)  y  K-|  (x) 
=  N[H-|  ](x),  obtenemos  el  siguiente  resultado  análogo  al  Teorema  3.5. 


Corolario  3.6. 


Sea  f  una  función  estrictamente  convexa  en  x0,  con  f'(x0)  *  0,  así  se  tiene 
(i)  Si  se  verifica  que 


n  +  1 


f"(x0) 

f'(x/ 


>  n 


Entonces  Hn  es  estrictamente  convexa  en  xQ  y  Hn  +  -|  es  estrictamente  cón¬ 
cava  en  xQ.  Por  tanto  f  es  n-  y  -convexa  y  no  es  (n+1)  -  y  -convexa. 


(¡i)  Si 


1  f"(x0) 

— 777 — J - 7  =  n 

v  (y o>  f'(xo)2 


,  entonces 


f  es  una  f  unción 


(n-1)  -  y  -convexa,  siendo  la  convexidad  de  Hn_i  y  la  concavidad  de  Hn+i 
estrictas  en  xQ.  # 


Por  otra  parte,  en  relación  con  la  curvatura,  se  prueba  que 
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Ejemplo  3.3. 

Consideremos  la  familia  de  funciones  convexas 

fm(x)  =  (1/2)  exp  (xm-1)  para  meR-{(0,1)}  y  x0  =  1.  Si  tomamos  cp(y)  =  log  y 
entonces  yQ  =  f(x0)  =1/2  y  vy(y)  =  (log  y)/2. 


En  estas  condiciones  obtenemos 


l- 


__1 _  f'm(xo) 


Al  variar  m  obtenemos  valores  diferentes  para  esta  cantidad.  Estos  valores 
nos  dicen  qué  términos  de  la  sucesión  {Hn}  son  funciones  convexas.  Por  ejemplo 
si  tomamos  m  =  -2,  (fig.1 9)  entonces  obtenemos  que  H1  y  H2  son  convexas  y  H3 
es  cóncava. 


HlW  =  2l'V"1'109  21 

log  2 

H2(X)  =  "  2 - 109  X 

H3(x)  =  ¿-|09  lj-  logx] 


K  (x)  = 


K2(X)  =  T~  '°gX 


Además  se  observa  (fig.  20.),  que 


K(f)(1)  >  K(H!  )(1)  >  K(H2)(1)  >  0  >  K(H3)(1). 
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Así,  podemos  definir 


Definición  3.7. 


Llamaremos  cp-grado  de  convexidad  de  f  en  x0  al  número  real  dado  por  la 
expresión 


(3.5  .) 


Vf](Xo>  =  I 


-1 

v"(y  o) 


f"(x0) 

f'(x/ 


si  xQ  no  es  un  mínimo  de  f .  En  otro  caso  U(p(f)(x0)  =  +  «>. 
Haciendo  unos  pequeños  cálculos  obtenemos 


(3.6)  lyfKxj 


<p'(yp)  f"(x  J 
<p"(y  r(x/ 


Notas  : 

Esta  cantidad  nos  da  la  medida  de  la  convexidad  antes  de  normalizar.  Por 
tanto  es  claro  que  la  normalización  anterior  no  afecta  nuestro  estudio,  solamente 
lo  simplifica. 

Además  podemos  notar  que  U^ÍXq)  e  [  0,  +  °°)  y  no  está  acotado,  como 
muestra  el  siguiente  ejemplo 


Ejemplo  3.4. 

Sea  la  familia  de  funciones  fn(x)  =(-1/n)  log  nx  en  [0,1],  con  ne  R+,  y  con¬ 
sideremos  la  familia  de  operadores  cpm(x)  =  2  +  x  -  mx2/2,  con  m  €  R+.  Si  toma¬ 
mos  xD  =1/n,  entonces  y0  =  f  (xQ)  =  0  y  U^Jf  ](^)  =  ^  . 
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De  la  Definición  3.7.  se  deducen  las  primeras  propiedades  de  este  índice  de 
medida  de  la  convexidad 


Corolario  3.8. 

Con  la  notación  anterior,  se  verifica  que 


(i)  Para  cada  número  real  r,  se  tiene 

\|/"(y  q) 

a)  Si  denotamos  k=— — - -  entonces 

V  (r+y  o) 

U9  [r+f](x0)  =  k  [f](x0). 

\|í"(y  ) 

b)  Si  r  >0  y  denotamos  k  = - ,  entonces 

V  (ryj 


U9  [r.f](x0)  =  (k/r)U(p[f](x0). 


(ii)  Si  1/f  es  convexa  en  (a,b),  entonces 

=  ~^Í7[2  +  v  "(yo>''ou/Kxo)  1  ■ 

y  O 

(iii)  Ujp  [f](x0)  >  0  .  Además,  [f](x)  =  0  en  W  sí  y  sólo  sí  f  es  afín  en  W,  es 
decir,  f(x)  =  mx  +  d  con  m,d  eR. 

(iv)  Sea  g  una  función  en  las  mismas  condiciones  que  f,  tal  que  gof  esté  defi¬ 
nida,  entonces 


Uv[f](xo)  . 
# 


i  V  (f(x0)) 

U  íg  of  ] ( x  )  =  U  [g](f  (x  ))  +  -  1  -  -  - 

•P  *  g(f(x0))  y  (g(f(x  ))) 
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3.3.  Grado  de  Convexidad. 

Finalmente  queremos  obtener  un  indice  de  convexidad  referente  a  todos  los 
operadores  cóncavos.  Para  ello  si  normalizamos  el  operador  9  por  la  transfor¬ 
mación 


<P*  (y)  = 


<p(y)  -  t<p'(y0)  +  <p"(y0)]  y 

[-  <P'(y0)] 


obtenemos  un  nuevo  operador  9*  de  las  mismas  características  geométricas  que 
9,  con  (  9*)’(y0)  =  1  y  (  9*)”(y0)  =-1  ■  Por  tanto,  de  (3.6.)  obtenemos 


Vf]<x0>  = 


rog 

f'(x0)2 


Así,  definimos  el  grado  de  convexidad  de  f  en  x0,  como  el  número  real  posi¬ 
tivo  dado  por  la  expresión 


r(Xo) 

(3.7.)  U[f](x  )  =  -~ 

f'(xn) 


si  x0  no  es  un  mínimo  para  f.  En  este  caso  U[f](x0)  =  +  °°. 

Es  fácil  probar  las  siguientes  propiedades  de  este  grado  de  convexidad 


Corolario  3.9. 

Se  verifica  que 

(i)  Para  todo  número  real  r  se  tiene 

a)  U[r  +  f](x0)  =  U[f]  (x0),  es  decir,  U[f]  es  invariante  por  traslaciones. 

b)  Si  r  e  R+,  entonces  U[r.f](x0)  =  U[f](x0)/r. 

(ii)  Si  1/f  es  convexa  en  (a,b),  entonces 
U[1/f](x0)  =  f(x0)  [2  -  f(x0)  U[f]  (x0)  ]. 
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(iii)  U[f]  (x0)  >  0  y  además  U[f](x)  =  0  en  (a,b)  sí  y  sólo  sí  f  es  afín,  es  decir, 
f(x)  =  mx  +  d  con  m,  d  eR. 

(iv)  Si  f  y  g  son  dos  funciones  crecientes,  se  verifica  que 
u[f  +  g](x0)  á  U[f]  (x0)  +  U[g]  (x0). 

(v)  Si  f  y  g  son  dos  funciones  crecientes,  tales  que  gof  está  definida,  entonces 

U[g  o  f](x0)  =  U[g]  (f(x0))  +  [1/g'(f(x0))]  U[f]  (xQ). 

Estos  resultados  nos  permiten  afirmar  que  U[f]  es  una  buena  medida  de  la  con¬ 
vexidad  en  cada  punto,  [19]. 


3.4.  Grado  de  Concavidad. 


Si  consideramos  una  función  cóncava  f  en  un  entorno  del  punto  x0,  y  cp  un 
operador  convexo  en  un  entorno  de  y0=f(x0),  podemos  realizar  un  proceso  simi¬ 
lar  al  que  hemos  hecho  para  definir  el  grado  de  convexidad,  y  obtenemos  así  el 
cp-grado  de  concavidad  de  f  en  xQ,  dado  por  la  expresión 


9'(yo) 

n 

9  9  (y0) 


[ 


'(xo) 


f  '(x/ 


] 


Como  indicamos  en  el  caso  del  cp-grado  de  convexidad,  normalizando 
rador  cp  obtenemos  un  índice  de  medida  de  la  concavidad  de  f  en  x0.  Este 
dado  por  el  número  real  positivo  obtenido  por  la  expresión 


el  ope- 
vendrá 


f  "(x0) 

(3.8.)  n  [f ](xQ)  = - 

f  (xo) 

si  x0  no  es  un  máximo.  En  este  caso  Pi  [f](x0)  =  +  °°. 

Geométricamente,  si  f  tiene  las  condiciones  indicadas  podríamos  definir  el 
grado  de  concavidad  de  f  en  xQ  como  el  grado  de  convexidad  de  -f  en  x0.  Ya  que 
de  (3.7.)  y  (3.8.)  obtenemos 

(3.9.)  n  [f ](X0)  =  U[-f](x0)  . 


Por  otra  parte,  se  sigue  de  (1 .9.)  que  el  grado  de  concavidad  de  f  en  x0  va  a 
tener  las  mismas  propiedades  que  el  grado  de  convexidad  de  f  en  xQ,  como 
hemos  visto  en  el  Corolario  3.9. 
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Uno  de  los  problemas  matemáticos  que  suelen  aparecer,  con  más  frecuen¬ 
cia  en  la  práctica,  es  la  resolución  de  una  ecuación  f(x)=0.  Aunque  en  algún  caso 
muy  sencillo  haya  procedimientos  de  resolución  directa  de  esta  ecuación,  en 
general,  nos  veremos  obligados  a  recurrir  a  procesos  iterativos  que  nos  permitan 
ir  obteniendo  valores,  que  se  espera,  sean  cada  vez  más  próximos  a  la  solución. 
A  partir  de  estos  valores,  mediante  nueva  aplicación  del  proceso  iterativo,  obte¬ 
nemos  otro  valor  más  próximo  a  la  solución  de  la  ecuación  y  así  sucesivamente. 
La  gran  mayoría  de  estos  procesos  iterativos  se  basan  en  la  transformación  de  la 
ecuación  original,  anteriormente  citada,  en  otra  equivalente  del  tipo  x=F(x).  La 
forma  de  obtener  F  depende  de  cada  método  particular.  Una  vez  obtenida  dicha 
función  se  construye  una  sucesión  de  la  forma  xn+-|=  F(xn),  la  cual  nos  debe  pro¬ 
porcionar  sucesivas  aproximaciones  a  la  solución  de  la  ecuación  f(x)  =  0  [7]  y  [12], 
Por  tanto,  transformamos  un  problema  de  resolución  de  una  ecuación  en  un  pro¬ 
blema  de  cálculo  de  puntos  fijos. 

En  el  primer  Capítulo  de  este  Tema  analizamos  los  métodos  iterativos  de 
Whittaker  y  de  Newton-Raphson  para  la  resolución  de  ecuaciones.  Comenzamos 
recordando  algunas  generalidades  acerca  de  dichos  métodos.  Posteriormente 
estudiamos  la  convergencia  de  estos  métodos  iterativos  utilizando  los  grados  de 
concavidad  exponencial  y  de  convexidad  logarítmico  respectivamente.  Además, 
mediante  estos  conceptos,  indicamos  posibles  procedimientos  para  aumentar  la 
velocidad  de  convergencia  de  dichos  métodos  iterativos. 

En  el  segundo  Capitulo  de  este  Tema  vamos  a  estudiar  la  convergencia  de 
algunos  procesos  iterativos  de  tercer  orden  [22],  es  decir  con  convergencia  cúbi¬ 
ca,  para  la  resolución  de  la  ecuación  no  lineal  f(x)  =  0.  Este  estudio  lo  realizamos 
a  partir  de  la  convexidad  de  la  curva  y  =  f(x).  Para  medir  esta  convexidad  utiliza¬ 
mos  el  grado  de  convexidad  logarítmico  de  la  función  f.  El  hecho  de  utilizar  esta 
medida  de  la  convexidad  se  debe  a  obtener  una  mayor  simplificación  operacional, 
si  bien  resulta  evidente  que  se  puede  realizar  un  estudio  análogo  utilizando  el 
grado  de  convexidad  de  f,  ya  que  Lf(x)  =  f(x)U[f](x). 
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El  resultado  de  Gander  [5]  en  el  que  se  prueba  que  todo  proceso  iterativo  de 
la  forma 


xn  +  1=  F(xn)  =  xn  -  H(Lf(xn))  f(xn)/f’(xn) 

donde  H(0)=1,  H’(0)  =  1/2  y  IH”(0)l  <  +«=,  tiene  convergencia  cúbica,  nos  lleva  a 
expresar  los  procesos  iterativos  de  tercer  orden  en  función  del  grado  de  convexi¬ 
dad  logarítmico  de  la  función  f.  Así,  se  puede  comprobar  que  los  procesos  iterati¬ 
vos  usuales  de  tercer  orden,  el  método  de  Chebyshev  [4]  y  el  método  de  Halley 
[5],  admiten  una  expresión  de  este  tipo.  Por  ello,  damos  resultados  de  conver¬ 
gencia  de  estos  métodos  imponiendo  condiciones  sobre  Lf  y  Lf>,  es  decir  sobre 
las  derivadas  segunda  y  tercera  de  la  función  f,  lo  que  es  habitual  para  asegurar 
la  convergencia  de  procesos  iterativos  de  tercer  orden. 
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CAPÍTULO  I 

LA  CONCAVIDAD  Y  LA  CONVEXIDAD  EN  LA  CONVERGENCIA 
DE  PROCESOS  ITERATIVOS  DE  ORDEN  MENOR  QUE  TRES 


1.1.  El  método  de  Whittaker 

1.1.1.  Generalidades 

El  método  de  Whittaker  es  un  proceso  iterativo  que  surge  como  consecuen¬ 
cia  del  Teorema  del  Punto  Fijo  [1 5].  Es  conocido  que  si  F:  [a,b]  [a,b]  es  una  apli¬ 

cación  contractiva,  y  partimos  de  un  punto  x0  arbitrario  de  [a,b],  la  sucesión  dada 
por  xn  =  F(xn .  i )  converge  al  único  punto  fijo  de  F.  Así,  si  consideramos  una  ecua¬ 
ción  de  la  forma  F(x)  =  x  equivalente  a  f(x)  =  0  en  [a,b]  trasladamos  el  problema 
de  resolución  de  la  ecuación  f(x)  =  0  a  un  problema  de  cálculo  del  punto  fijo  de  un 
operador  F.  Es  claro  que,  si  f:  [a,b]  R  es  una  aplicación  que  se  anula  para  un 
único  punto  de  [a,b],  podemos  definir  F:  [a,b]  R  de  manera  que  F(x)  =  x-Xf(x), 
para  X  *  0  y  la  solución  de  f(x)  =  0  será  el  único  punto  fijo  del  operador  F.  Así  el 
método  de  Whittaker  consiste  en  determinar  valores  de  X,  distintos  de  cero,  para 
los  cuales  la  aplicación  F  sea  contractiva  y  F(x)  pertenezca  al  intervalo  [a,b]  para 
todo  punto  x  e  [a,bj.  Entonces  estamos  en  disposición  de  hallar  el  único  punto  fijo 
de  F  en  [a,b],  que  será  la  única  raíz  de  f  en  dicho  intervalo. 

Así,  como  hemos  indicado  anteriormente,  construimos  la  sucesión  {xn}  dada 
por  la  expresión 

(1.1)  xn  =  F(xn  _  i )  =  xn  _  -  A.f(  xn  _  i )  conx0e[a,bj. 

A  continuación  nos  planteamos  el  estudio  de  condiciones  suficientes  para 
que  la  sucesión  {xn},  dada  por  (1.1),  sea  convergente  a  s,  única  raíz  de  la  ecua¬ 
ción  f(x)  =  0  en  [a,bj. 

El  resultado  siguiente  nos  proporciona  condiciones,  que  debe  verificar  f,  para 
poder  aplicar  el  método  de  Whittaker,  obteniendo  una  sucesión  {xn}  convergente 
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a  s.  Además  nos  indica  cómo  determinar  los  valores  de  k  para  los  cuales  se  veri¬ 
fica  dicha  convergencia,  su  demostración  puede  verse  en  [11], 

Teorema  1.1. 

Sean  f  e  c(  ^([a,b]),  f(a)f(b)  <  0,  f'(x)  *  0  en  [a,b],  {xn}  la  sucesión  dada  por 
(1.1),  con  f(x0)  >  0,  M  =  máx.{lf'(x)l  |  x  e[a,b]}  entonces 

i)  Si  f  es  creciente  tenemos  que 

a)  para  k  e  (  0,  1/M  ]  la  sucesión  {xn}  es  decreciente  a  s. 

b)  para  ke  (  1/M,  2/M  )  la  sucesión  {xn}  converge  a  s  pero  no  sabemos  su 
forma. 

ii)  Si  f  es  decreciente  tenemos  que 

a)  para  ke  ( -1/M,  0  )  la  sucesión  {xn}  es  creciente  a  s. 

b)  para  k  e  ( -2/M,  -1/M  ]  la  sucesión  {xn}  converge  a  s  pero  no  sabemos 
su  forma.  # 

Notar  que  la  elección  de  x0e  [a,b]  puede  hacer  variar  el  carácter  creciente  o 
decreciente  de  la  sucesión. 

Otro  tipo  de  condiciones  de  convergencia,  así  como  el  análisis  del  error  y 
orden  de  convergencia  del  método  pueden  verse  en  [13]. 

El  Teorema  1.1  nos  proporciona  los  valores  de  k  para  los  que  la  sucesión  {xn} 
dada  por  (1.1)  es  convergente.  A  continuación,  queremos  determinar,  en  los 
casos  en  que  la  sucesión  sea  decreciente  o  creciente,  los  valores  de  k  que  nos 
proporcionen  una  situación  óptima  del  método  de  Whittaker.  Así,  siguiendo  con  la 
notación  de  dicho  teorema,  obtenemos  el  siguiente  resultado 

Teorema  1 .2. 

Sean  k  y  p  dos  números  reales,  xn  =  xn_i  -  k  f(xn  _  -j )  =  F(xn_-|)  e  yn= 

yn-i  -ftf(yn-i)  =  G(yn-i)con  xo  =  yo 

i)  Si  f  es  creciente  en  [a,b]  y0<p<A.<1/M  entonces  la  sucesión  {xn  }  con¬ 
verge  a  s  más  rápidamente  que  la  {yn}. 

ii)  Si  f  es  decreciente  en  [a,b]  y  -1/M  <  k  <  p  <  0  entonces  la  sucesión  {xn} 

converge  a  s  más  rápidamente  que  la  {yn}.  # 

La  demostración  de  este  resultado  quede  verse  en  [11],  Como  consecuencia, 
podemos  deducir  que  a  mayores  valores  de  IA.I  la  sucesión  obtenida  tiene  mayor 
velocidad  de  convergencia. 
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1. 1.2.  La  Concavidad  en  la  Convergencia  del  método  de  Whittaker 

A  continuación  vamos  a  efectuar  un  estudio  del  método  de  Whittaker  median¬ 
te  el  grado  de  concavidad  exponencial  de  f.  Para  ello  consideramos  en  todo 
momento  fe  c(  m([a,b])  con  m  >  2,  f  una  función  decreciente  y  cóncava.  Si  f  fuera 
creciente  bastaría  con  efectuar  el  cambio  de  f(x)  por  f(-x). 

También  elegimos  x0  e  [a,b]  con  f(x0)  >  0,  en  otro  caso  los  resultados  serán 
análogos. 

Comenzamos  este  estudio  probando  un  resultado  de  convergencia  global 
para  el  proceso  iterativo  de  Whittaker,  en  función  del  grado  de  concavidad  expo¬ 
nencial  de  f.  Si  bien  este  resultado  no  posee  gran  interés  práctico,  dado  que  sus 
hipótesis  son  excesivas  respecto  del  Teorema  1 .1 ,  sí  posee  interés  teórico  ya  que 
dichas  hipótesis  se  expresan  mediante  el  grado  de  concavidad  exponencial  de  la 
función.  Además  este  resultado  acota  los  posibles  valores  del  parámetro  X  en  fun¬ 
ción  del  grado  de  concavidad  exponencial  de  la  función  f. 


Teorema  1.3. 

Sea  f  en  las  condiciones  indicadas  verificando  que  1/lf’(a)l<b-a,  si  X  e  R_  y 
Ef(x)  >  M'(b)  para  x  e  [a,b]  entonces  la  sucesión  {xn}  dada  por  (1.1),  con 
x0  e  [a,b]  y  f(xo)>0,  converge  de  forma  creciente  a  la  única  raíz  de  f(x)  =  0  en  [a,b]. 

Demostración: 

La  estructura  de  esta  demostración  va  a  ser  análoga  a  la  del  Teorema  1.1. 
Aquí  f(a)  >  0  >  f(b). 

a)  Probemos,  por  inducción  sobre  n,  que  xn  <  s  para  cualquier  neN.  Como 
f(x0)  >  0,  entonces  x0  <  s  y  por  tanto  x-|  -  s  =  F(  x0)  -  F(s),  con  F(x)  =  x  -  Á.f(x). 

Por  el  Teorema  del  Valor  medio,  existe  0oe  (s,  x0)  tal  que  x-j  -  s  =  F’(  0O)  (  x0-  s) 
=[1-  A.f’(  0O)](  x0-  s  ). 

Si  probamos  que  A.f’(  0O)  <  1 ,  como  x0-  s  <  0,  es  claro  que  x-|  -  s  >  0.  Para 
probar  esto,  como  Ef  (x)  >  M’(b)  para  x  e  [a,b],  se  tiene  que 

b  b 

1  -•M*"  •  vk  ~  ímm  =l,,(b)  <b‘a) 

3  3 
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por  otra  parte 


f'(b)  f'(a) 


1 

lf'(a)l 


1 


1 


lf'(b)l  lf'(a)l 


<  b  -  a, 


pues  0  >f'(a)  >f'(b), 


luego  X  f’(b)(b-a)  <  b-a. 

Entonces  como  f’  es  decreciente,  f’(b)  <  f’(0o)  y  X  <  0  se  sigue  que 


Xf'(0o)<W'(b)<^J  =  1 

Suponiendo  cierto  que  xn_  ^  <  s,  probémoslo  para  xn.  xn-s  =  F(xn.  i)  -  F(s) 
=  F’(en .  1  )(xn .  1  -  s)  =  [1  -  W’(en  - 1  )Kxn  - 1  -  s)  para  en  - 1  €  (s’xn  - 1  )■  Como  en  el 
caso  anterior,  basta  probar  que  M’(  6n  - 1 )  <  1  y  como  Xf’(  0n .  )  <  Xf’(b)  <  1 ,  se 
sigue  que  xn<  s. 

b)  Veamos  que  la  sucesión  {xn}  es  decreciente 

xn .  i  -xn  =  xn  - 1 )  <  0  al  ser  ^  <  0  y  f(  xn  - 1 )  >  °- 

c)  Luego  la  sucesión  {xn}  es  creciente  y  acotada,  por  tanto  tiene  limite,  sea 

este  p.  Entonces  pasando  al  límite  en  (1.1)  tenemos  que  p  =  p  -  A.f(p),  luego  f(p) 
=  0  y  por  tanto  p=s.  # 

Como  hemos  comentado  anteriormente,  este  teorema  sólo  tiene  un  interés 
teórico  pues  la  acotación  obtenida  X  >  Ef  (x)/f ’(b)  es  peor  que  la  obtenida  ante¬ 
riormente  si  Ef  (x)  <  1 .  Por  otra  parte,  en  el  caso  de  que  Ef  (x)  >1  va  a  ser  incom¬ 
patible  con  las  hipótesis  del  teorema,  ya  que  si  consideramos  f  como 
habitualmente  y  suponemos  que  Ef(x)  =  k  >  1,  entonces 


b-  a  < 


b 


k(b-  a)  =  J*Ef (t  )dt 

a 


1 

p(b) 


i  i 

f'(a)  <lf'(a)l  ’ 


de  lo  que  podemos  deducir  que  no  se  puede  encontrar  un  intervalo  [a,b]  en  el  que 
1/lf’(a)l  <  b-a  y  a  la  vez  el  grado  de  concavidad  exponencial  de  f  sea  mayor  que 
uno.  A  continuación  damos  un  ejemplo  de  esta  situación,  así  como  del  caso  en 
que  sí  podemos  aplicar  el  teorema,  poniendo  de  manifiesto  que  obtenemos  peor 
acotación  para  los  valores  del  parámetro  X  que  en  el  Teorema  1.1. 
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Ejemplo  1.1. 


(i)  Estudiemos  la  función  f  dada  por  f(x) 
f  e  c(2((-°°,1/2))  y  [a,b]  1  1/2)  tal  que  a<0  y  t»0 


log(1-2x)/2,  (fig.21) 


f‘(x) 


1  . 
2x  -  1  ’ 


f"(x)  = - E  (x) 

(2x  -  D 


f"(x)  _2 
(f'(x))2 


Vemos  que  la  función  es  cóncava,  decreciente  y  que  su  grado  de  concavidad 
es  dos.  Por  tanto,  según  lo  visto  en  la  nota  previa,  no  se  podrá  aplicar  el  teorema 
anterior.  En  este  caso,  el  que  1/lf’(a)l<b-a  es  equivalente  a  que  1  <  b  +  a  y  como 
a  <  0  se  sigue  que  b  >  1 ,  con  lo  cual  no  es  posible,  como  ya  habíamos  indicado, 
aplicar  el  teorema  anterior. 


(ii)  Consideremos  la  función  f(x)  =  e2  -  ex  y  [a,b]  =  [1,3].  (fig.22.). 

Entonces  es  claro  que  1/lf’(a)l  =  e‘a  <  2  =  b  -  a  y  Ef(x)  =  1/ex  <  1  en  [1,3], 
luego  si  M’(b)  <  Ef(x)  en  [1 ,3],  es  decir,  X{-  e2)  <  e"x  para  x  e  [1 ,3].  Entonces,  por 
el  resultado  anterior,  el  método  de  Whittaker  es  convergente  a  la  raíz  s  =  2.  Ahora 
bien,  según  lo  indicado  ha  de  verificarse  que  0  >  X  >  -e'2,  mientras  que  el  Teorema 
1.1.  nos  permite  considerar  -2  e"  2<  X  <0.  Luego  la  acotación  dada  por  el  Teorema 
1 .3  es  menos  fina  que  la  que  nos  proporciona  el  Teorema  1 .1 ,  teniendo  en  cuen¬ 
ta  que  el  método  de  Whittaker  posee  mayor  velocidad  de  convergencia  para 
mayores  valores  de  IA.I. 
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La  convergencia  del  proceso  iterativo  de  Whittaker  es  lineal.  A  continuación 
estudiamos  la  influencia  que  tiene  la  concavidad  de  la  función  f  en  la  velocidad  de 
convergencia  de  la  sucesión  {xn},  estudio  que  realizamos  mediante  el  grado  de 
concavidad  exponencial  de  la  función. 


Teorema  1.4. 


Si  g  es  una  función  en  las  mismas  condiciones  que  f,  siendo  Eg(x)  <  Ef  (x)  en 
[a,b],  entonces  existe  k  =  máx.{1,f’(b)/g’(b)}  de  manera  que  la  función  h(x)  =  kg(x), 
nos  proporciona  una  sucesión  yn  =  yn .  1  -X  h(yn .  i ),  con  y0  =  x0,  que  converge 
de  forma  creciente  a  s  para  0  >  X  >  máx.{  [  g'(b)  ]_1 ,  [f'(b)  ]  }.  Además  se  veri¬ 

fica  que  {yn}  converge  más  rápidamente  a  s  que  {xn},  es  decir  yn  >  xn  para  cual¬ 
quier  neN. 

Demostración  : 


Como  f  y  g  son  cóncavas  y  decrecientes,  las  sucesiones  {xn}  e  {yn}  conver¬ 
gen  a  s  y  son  crecientes  para  el  valor  de  X  considerado.  Entonces  veamos  que 
yn  >  xn  por  inducción. 


Supongamos  que  g'(b)  <  f'(b),  x-|  -  y-|  =X  [g(x0)  -  f(x0)],  ahora  bien  Eg(x)  < 
Ef(x)  en  [a,b]  luego  se  verifica 


J[Ef(t )  -  Eg(t  )]dt  >  0  «. 

X 


f'(x)  >  g'(b) 


1 

f'(b) 


>0. 
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Entonces  g'(x)  <  f'(x)  en  (a,b),  de  donde  se  sigue  que  para  x  e(a,s) 
s 

i  [f '  (t )  —  g‘  (t )  ]dt  >0  ,  luego  g(x)  >  f  (x)  en  (a,s) 
x 

y  por  tanto  g(x0)  >  f(x0)  y  como  X  <  0  se  sigue  que  x-)  <  y-j . 

Supongamos  que  y^  >  x^  para  k=1 ,2,...,  n-1 ,  entonces  si  llamamos  G(x)  =  x  - 
X  h(x),  como  F’(x)  >  0  al  ser  X  >  1/f ’(b)  tenemos  que  xn-  yn  =  F(xn .  i )  -  G(yn .  i )  < 
F(yn  - 1 )  -  G(yn .  -i )  =  X[g(yn .  i )  -  f(xn  _-,)]<  0,  al  ser  X<  0  y  g(yn .  >  f(xn .  luego, 

xn  <  yn- 

Si  g’(b)  >  f’(b),  como  k  =  f’(b)  [g’(b)]‘  ,  consideraremos  h(x)  =  k  g(x)  que  veri¬ 
fica  Eh(x)  =  g’(b)  [f’(b)]‘  ■’Egíx)  <  Eg(x)<Ef(x)  en  [a,b]  por  ser  g'(b)  [f*(b)]“ 1  <  1. 
Además  h’(b)  =  f’(b),  entonces  aplicando  la  situación  anterior  a  la  función  h  queda 
probado  el  resultado.  # 


Ejemplo  1.2. 

Consideramos  las  funciones  f(x)  =  1-  x3/216,  g(x)=1-  x2/36.  (fig.  23.).  Vemos 
que  ambas  tienen  las  mismas  características  en  [5,7]  ,  además  de  la  misma  raíz 
s=6.  Tenemos  que 

f’(x)  =  -  x2/72;  f”(x)  =  -  x  /36;  Ef  (x)  =  144/x3 

g’(x)  =  -  x/18;  g”(x)  =-1/18;  Eg(x)  =  18/x2 

luego  las  dos  funciones  son  cóncavas  y  decrecientes  en  dicho  intervalo  ,  además 
se  verifica 

Eg(x)  =  1 8/x2  <  1 44/x3  =  Ef  (x)  <=>  1 8  x  <  1 44  o  x  <  8 

Por  otra  parte  f’(b)  =  f’(7)  =  -  49/72  y  g’(b)  =  g’(7)  =  -7/18  luego  k  =  7/4  >  1 . 
Vemos  que,  aún  teniendo  menor  grado  de  concavidad  la  función  g  que  f,  la  apro¬ 
ximación  es  mejor  con  f.  (Tabla. 1.) 

Esto  se  corrige,  procediendo  como  indica  el  teorema  anterior,  definiendo  la 
función  h(x)  =  f’(7)  [g’(7)]‘  1  g(x)  =  (7/4)  g(x).  (Tabla.  2. ). 

También  es  interesante  notar  que  la  relación  dada  entre  los  grados  de  con¬ 
cavidad  exponencial  de  f  y  g,  sólo  es  necesaria  en  (a,s),  es  decir,  donde  se  toma 
x0.  Ahora  bien,  como  s  es  un  punto  desconocido,  lo  exigiremos  en  todo  [a,b], 
teniendo  en  cuenta  la  posible  variabilidad  del  punto  b 
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TABLA  .  1  . 


«nai 

f 

9 

i 

7 

1 

2 

6.136054421768700 

6.4693877551021 

3 

6.033812309317900 

6.2304906969036 

4 

6.008830357183800 

6.1 154288877524 

5 

6.002333195155500 

6.05834845995 

6 

6.000618344285500 

6.0296306614936 

7 

6.000164003704800 

6.0150818486762 

8 

6.000043507893500 

6.0076855365834 

9 

6.00001 1542678700 

6.0039187812172 

na 

6.000003062326900 

6.001998751362 

WTU 

6.000000812453000 

6.0010196080416 

wvm 

6.000000215548600 

6.0005201657518 

wrm 

6.000000057186400 

6.000265379646 

wrm 

6.000000015171900 

6.000135394904 

ra 

6.000000004025100 

6.0000690782844 

V _ J 
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TABLA . 2  . 


E39 

h 

i 

7 

2 

6.0714285714286 

3 

6.0098396501456 

4 

6.0013987486841 

5 

6.0001996814908 

6 

6.0000285230792 

7 

6.0000040746674 

8 

6.0000005820942 

9 

6.0000000831563 

MÜM 

6.00000001 18794 

m 

6.0000000016971 

m 

6.0000000002424 

wrm 

6.0000000000346 

wrm 

¿  bóbóbódóoóbé 

ra 

6.0000000000006 

Otra  cuestión  interesante,  como  ya  vimos  anteriormente,  es  la  optimización 
del  método.  Es  conocido,  que  a  mayor  valor  del  parámetro  en  valor  absoluto,  el 
método  proporciona  una  más  rápida  aproximación.  Este  resultado  lo  podemos 
obtener  como  una  consecuencia  inmediata  del  teorema  anterior. 


Corolario  1.5. 

Sean  X  y  p  dos  números  reales  tales  que  [f ’(b)]"  1  <  X  <  p  <  0  si  xn  =  xn .  i  - 
M(xn .  i )  e  yn  =  yn  _  i  -  n  f(yn .  1 )  con  xQ  =  y0  y  f(xQ)  >  0,  entonces  yn  <  xn  para 
cualquier  n  e  N. 

Demostración  : 

Basta  considerar  la  función  g(x)  =  (|x/^.)f(x),  la  cual  verifica  las  condiciones  del 
teorema  anterior  y  como  g’(b)  <  f’(b),  entonces  tenemos  que  {xn}  converge  a  s 
más  rápidamente  que  {yn},  es  decir,  yn  <  xn  para  cualquier  n  e  N.  # 
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1 .2.  El  método  de  Newton-Raphson. 

1.2. 1.  Generalidades 

La  idea  de  construcción  del  método  consiste  en  encontrar  una  aproximación 
de  una  raíz  s  de  f,  resolviendo  para  ello  la  ecuación  lineal  L(x)  =  0,  donde  L  es  el 
polinomio  de  interpolación,  en  este  caso  de  Taylor  [18],  construido  a  partir  de  un 
punto  inicial  xQ  el  valor  de  la  función  y  el  valor  de  la  derivada  en  dicho  punto.  Así 
el  polinomio  L(x)  =  f(x0)  +  (x-x0)f’(x0)  interpola  a  f  en  un  entorno  de  x0,  ahora 
igualando  a  cero  y  despejando  tenemos 


f(x0) 

x  =  x_ - . 

0  f'(x0) 

Notemos  que  ha  de  verificarse  que  f’(x0)  *  0.  Tomando  xi  =  x  y  repitiendo  el 
proceso  obtenemos  una  nueva  aproximación  de  s  que  denotaremos  X2,  de  esta 
manera  conseguimos  una  sucesión  {xn},  definida  por  el  proceso  iterativo 

<12» 

Newton  fue  el  que  originalmente  dio  una  idea  acerca  de  este  método,  pero 
Raphson  expresó  la  fórmula  de  recurrencia  por  vez  primera,  por  lo  cual  a  este 
método  se  le  denomina  de  Newton-Raphson. 

La  construcción  de  xn  +  i  a  partir  de  xn  se  puede  ver  como  el  punto  de  corte 
de  la  tangente  a  la  curva  y  =  f(x)  en  el  punto  (xn,f(xn))  con  el  eje  OX.  Es  por  ello 
que  también  se  conoce  este  proceso  iterativo  como  el  método  de  la  tangente. 

De  momento  hemos  construido  una  sucesión  {xn},  pero  no  sabemos  si  con¬ 
verge  a  la  raíz  s  de  f(x)  =  0.  Daremos  a  continuación  algunos  resultados  que  nos 
aseguren  la  convergencia  del  método.  La  demostración  del  primero  de  ellos 
puede  verse  en  [16].  Denotaremos  por  <x,y>  el  intervalo  más  pequeño  que  con¬ 
tiene  a  ambos  puntos,  es  decir,  [x,y]  ó  [y,x]. 
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Teorema  1 .6. 

Sea  f  :  R-»R.  Consideramos  h0=  -f(x0)/f’(x0)  y  un  intervalo  J0  =  <x0,x0+ 
2h0>,  cumpliendo  que  f(x0)f’(x0)  *  0  y  en  el  que  exista  f”.  Sea  también  M  =  sup. 
{ |  f”(x)  |  I  x  e  J0  }.  Entonces  si 

(1.3)  2 1  h0 1 M  <  |  f’(x0)  l 

se  cumple  que,  tomando  como  punto  inicial  x0,  la  sucesión  {xn}  definida  por 
(1 .2)  está  contenida  en  JQ  y  el  límite  de  la  misma,  cuando  n  tiende  a  infinito,  es  s 
(  único  cero  en  J0  de  f ).  Además,  si  s  *  x0+  2h0  _  se  tienen  las  siguientes  rela¬ 
ciones  : 


(a) 

(b) 


I  x  -  xn  I 

n+1  n 


M 


I  x  -  x  .  I 

n  n- 1 


2  2lf'(xn)l 


'S  '  Xn+1 '  -  2  If '  (xn)  I 


para  n  >  1 


2 

x  I  para 

n- 1  r 


n  >  1 . 


# 

Este  teorema  prueba  la  convergencia  del  método  sin  tener  en  cuenta  el  cono¬ 
cimiento  previo  de  la  raíz  s.  En  cambio,  el  siguiente  teorema  exige  el  conocimien¬ 
to  de  la  existencia  de  una  solución  s  de  f(x)  =  0. 

Teorema  1.7. 

Sea  f  tal  que  f”  sea  continua  y  f’(x)  *  0  en  J,  intervalo  abierto  que  contiene  a 
una  raíz  s  de  f.  Entonces  existe  e  >  0  cumpliendo  que  para  todo  punto  inicial,  per¬ 
teneciente  al  entorno  de  centro  s  y  radio  e,  el  método  de  Newton-Raphson  es  con¬ 
vergente. 

Demostración  : 

Sea  G(x)  =  x-f(x)/f’(x),  como  G’(s)  =  0  <  1,  como  consecuencia  del  Teorema 
del  Punto  Fijo  [15]  obtenemos  que  existe  e  >  0  tal  que  la  sucesión  xn+  i=  G(xn) 
converge  al  punto  fijo  de  G,  que  es  s.  # 


Teorema  1 .8. 

Sea  f(x)  :  [a,b]  ->  R  cumpliendo: 
(a)  f(a)f(b)  <  0 
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(b)  f’(x)  íOVxe  [a,b] 

(c)  f”(x)  no  cambia  de  signo  en  [a,b] 

entonces  existe  una  única  raíz  s  de  f(x)  =  0  en  [a,b],  a  la  cual  converge  la  suce¬ 
sión  {xn},  definida  por  xn+  -j -  xn-  f(xn)/f’(xn),  para  x0  e  [a,b]  tal  que  f(x0)f”(x0)  >  0. 

Para  su  demostración  ver  [7],  # 


Notas: 

1.  Las  hipótesis  de  este  teorema  reciben  el  nombre  de  condiciones  de 
Fourier,  bajo  las  cuales  podemos  asegurar  la  convergencia  del  método  de 
Newton-Raphson. 

2.  Si  imponemos  la  condición  |f(c)  |/|f’(c)  I  <  b  -  a,  donde  c  es  el  extremo 
de  [a,b]  tal  que  |f’(c)  |  es  mínimo,  entonces  se  puede  tomar  cualquier  punto  en 
[a,b]  como  punto  inicial.  En  nuestro  caso  c  =  a,  entonces  de  |f(c)  |/|f’(c)  I  <  b- 
a  se  tiene  que  b  -  a  >  lf(a)l/lf(b)l,  y  como  f(a)  <  0  y  f’(a)  >  0,  se  sigue  que  b  -  a 
>  (-f(a))/f’(a).  Ahora  si  tomamos  como  punto  inicial  x0  e  [a,s],  se  cumple  que  -f(x0) 
<  -f(a)  y  f  (xQ)  >  f’(a),  entonces  b  -  a  >  (-f(a))/f’(a)  >  (— f(x0))/f’(x0)  =  x-|  -  x0  =>  b  > 
x-|  -  (xQ-  a)  >  x-| ,  es  decir  que  x-|  e  [a,b]. 

Para  terminar  las  generalizaciones  relativas  al  método  de  Newton-Raphson, 
indicar  que  es  un  proceso  iterativo  cuyo  orden  de  convergencia  es  al  menos  cua- 
drático  [12], 

Un  estudio  más  detallado  de  este  proceso  iterativo  puede  verse  en  [8], 

1.2.2.  La  Convexidad  en  el  Método  de  Newton-Raphson 

A  continuación  vamos  a  efectuar  un  estudio  del  Método  de  Newton-Raphson 
en  función  del  grado  de  convexidad  logarítmico  de  f.  En  este  párrafo  considerare¬ 
mos  f  una  función  creciente,  en  el  caso  de  ser  una  función  decreciente  el  cambio 
de  f(x)  por  f(-x)  nos  conduce  a  la  situación  considerada.  Además  tomaremos 
x0€[a,b]  con  f(x0)  >  0. 

A  partir  de  la  interpretación  geométrica  del  método  de  Newton-Raphson, 
aparece  la  “necesidad”  de  que  f”(x)  no  cambie  de  signo  en  el  intervalo  [a,b]  para 
asegurar  la  convergencia.  Entonces  nos  preguntamos  qué  ocurre  si  f”  cambia  de 
signo  en  [a,b].  En  principio,  pueden  presentarse  dos  situaciones.  En  primer  lugar, 
que  la  raíz  de  f(x)  =  0  sea  un  punto  de  inflexión  de  f  en  [a,b],  entonces  tanto  si  el 
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punto  inicial  x0  está  a  la  derecha  como  si  está  a  la  izquierda  de  la  raíz,  habrá  que 
asegurar  que  la  sucesión  {xn}  esté  contenida  en  [a,b],  En  segundo  lugar,  que  la 
raíz  de  f(x)  =  0  no  sea  un  punto  de  inflexión  de  f  en  [a,b].  Aquí,  puede  ocurrir  que 
el  punto  inicial  x0  esté  a  la  derecha  o  a  la  izquierda  del  punto  de  inflexión.  Si  por 
ejemplo  el  punto  de  inflexión  está  a  la  derecha  de  la  raíz,  el  punto  inicial  x0  puede 
estar  comprendido  entre  la  raíz  y  el  punto  de  inflexión,  en  cuyo  caso  la  conver¬ 
gencia  de  la  sucesión  {xn}  a  la  raíz  está  asegurada  por  el  Teorema  1 .8.  de  con¬ 
vergencia.  Si  por  el  contrario,  el  punto  inicial  x0  está  a  la  derecha  del  punto  de 
inflexión,  entonces  habrá  que  asegurar  que  la  sucesión  {xn}  esté  contenida  en 
[a,b].  Situaciones  análogas  se  pueden  presentar  si  el  punto  de  inflexión  se 
encuentra  a  la  izquierda  de  la  raíz. 

Todas  estas  observaciones  muestran  que  las  condiciones  de  Fourier  son 
insuficientes  para  asegurar  la  convergencia  del  método  de  Newton-Raphson  en  el 
caso  de  que  f”  cambie  de  signo,  y  por  tanto,  la  función  tenga  un  punto  de  inflexión 
en  su  dominio.  A  continuación  obtenemos  un  nuevo  resultado  de  convergencia 
para  el  método  de  Newton-Raphson  que  generaliza  los  resultados  de  convergen¬ 
cia  conocidos. 


Lema  1.9. 

Si  {xn}  es  la  sucesión  dada  por  (1.2),  x0e[a,b],  ILf(x)l  <  1  en  [a,b],  con  F(x) 
dada  por  (1.2),  y  a  <  F(x0)  <  b,  entonces  {xn}  c  [a,b]. 

Demostración: 

Obviamente  x1  =F(x0)e[a,b],  Además,  por  el  Teorema  del  Valor  Medio,  tene¬ 
mos  que  x2-s  =  F(x-|  )-F(s)  =  F’(x-)  )(x-|  -s)  con  x1  e(s,x1 ).  Por  otra  parte,  como 
IF’(x)l<1  en  [a,b]  entonces  existe  0<M<1  tal  que  IF’(x)l<M  en  [a,b].  Así,  lx2-sl  < 
M.lx1  -si  y  aplicando  inducción  obtenemos  que  lxn-sl  <  Mn"  ^  Ix-j  -si  <  lx-(  -si.  Por 
lo  tanto,  xne  [a,b]  para  todo  ne  N.  # 

Teorema  1.10. 

En  las  condiciones  del  Lema  anterior,  si  feC(m([a,b])  con  m>2,  f(a)f(b)<0  y 
f’(x)*0  en  [a,b],  entonces  el  proceso  iterativo  dado  por  (1.2)  es  convergente  para 
cualquier  x0e[a,b]. 

Demostración: 

Como  en  el  lema  previo,  sabemos  que  existe  0  <  M  <  1  tal  que  lxn-sl<Mnlxo-sl 
en  [a,b]  para  todo  neN,  luego  lím  lxn-sl=0.  Por  lo  tanto,  {xn}  es  una  sucesión  con¬ 
vergente  as.  # 


81 


LA  CONVEXIDAD  EN  LA  RESOLUCIÓN  DE  ECUACIONES  ESCALARES  NO  LINEALES 


Ejemplos  1.3. 

1)  Dada  la  función  f(x)=-x3+3x2-2,  buscamos  una  aproximación  a  la  solución 
de  la  ecuación  escalar  f(x)=0  en  el  intervalo  [j^  ,  -j^].  Esta  puede  obtenerse 
mediante  la  sucesión  numérica  {xn}  dada  por  el  método  de  Newton. 

2  - 

1  • 

0.25  0.5  0.75  /1  1.25  1.5  1.75 

Fig.  24. 

-2¥  — 

(6(-x3+3x2-2)(1-x) 

Es  fácil  obtener  que  F’(x)  =  2  R  \  y  Por  tanto  comprobar  que 

ri  i9i  (-vjx  +bx) 

ILf(x)l<1  en  L]~q  ,  Para  comprobar  que  la  sucesión  {xn}  está  contenida  en 

[1  191  1  19 

ÍO  ’  To-I’  es  suficiente  ver  por  el  Lemal.9  que  ^x^y^.  Tomando  por  ejemplo 

como  punto  de  salida  xg=1’6,  es  claro  que  x-|=0775e[^Q  ,  -j^].  Entonces  por  el 

Teorema  1.10  resulta  que  el  método  de  Newton  nos  proporciona  una  sucesión 

convergente  a  la  raíz  de  la  ecuación  escalar,  como  puede  verse  en  la  Tabla  3. 


iter. 

xn 

0 

1 .6000000000000000 

1 

0.7750000000000000 

2 

1.0079986833443050 

3 

0.9999996588133421 

4 

1 .0000000000000000 

Tabla  3 


Por  otra  parte,  ver  Fig.  24.,  la  función  presenta  un  punto  de  inflexión  en  la 
raíz.  Por  ello,  esta  función  no  verifica  las  condiciones  de  Fourier  y  en  consecuen¬ 
cia  no  podemos  asegurar  su  convergencia  con  las  hipótesis  habituales  de  con¬ 
vergencia  [16]  para  este  método. 
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2)  Dada  la  fundón  f(x)  =  ^  +  sen  x,  buscamos  una  aproximación  a  la  solución 
de  la  ecuación  no  lineal  f(x)=0  en  el  intervalo^  ,  mediante  el  método  de 
Newton.  Se  comprueba  fácilmente  que  ILf(x)l<1  en  [-1  '00297, 0’634867]c[:^  ,  |], 
Por  lo  tanto,  restringimos  el  intervalo  inicial  al  intervalo  [-1  '00297, 0’634867],  en  el 
cual  se  verifican  todas  las  condiciones  del  Teorema  1.10.  Tomando  por  ejemplo 
como  punto  de  salida  xo=0’6,  es  claro  que  x-|=-0’6899e[-1  '00297, 0’634867], 


Por  lo  tanto,  por  el  Lema  1 .9,  se  sigue  que  la  sucesión  {xn}  está  contenida  en 
[-1’00297,0’634867),  y,  por  el  Teorema  1.10,  el  método  de  Newton  nos  proporcio¬ 
na  una  sucesión  convergente  a  la  raíz. 


iter. 

xn 

0 

0.6000000000000000000 

1 

-0.6899509655978506667 

2 

-0.51 297262471 50719697 

3 

-0.5235667752006047706 

4 

-0.5235987753027045709 

5 

-0.5235987755982988737 

6 

-0.5235987755982988705 

7 

-0.5235987755982988742 

8 

-0.523598775598298871 1 

9 

-0.523598775598298871 1 

Tabla  4 


Notemos  sin  embargo,  Fig.  25.,  que  para  poder  aplicar  las  condiciones  habi¬ 
tuales  de  convergencia,  tendríamos  que  restringirnos  al  intervalo  ,  a]  con  a<0. 


83 


LA  CONVEXIDAD  EN  LA  RESOLUCIÓN  DE  ECUACIONES  ESCALARES  NO  LINEALES 


Por  tanto,  hemos  aumentado  el  dominio  de  los  puntos  de  salida  para  la  aplicación 
del  método  de  Newton. 

En  las  condiciones  de  Fourier,  para  f”(x)  >  0  y  x0e  [a,b]  con  f(x0)  >  0,  por  el 
Teorema  1 .10  se  prueba  que  la  sucesión  {xn},  dada  por  (1 .2),  es  decreciente  a  la 
raíz  s.  En  esta  situación  estudiamos  a  continuación  la  influencia  que  tiene  la  con¬ 
vexidad  de  la  función  en  la  velocidad  de  convergencia  del  método  iterativo. 

Sea  g  una  función,  en  las  mismas  condiciones  que  f  en  [a,b],  tal  que  g(s)  =  0. 
Consideramos  la  sucesión  yn  =  G(yn.-j),  siendo  G(y)  =  y  -  g(y)/g’(y)  ey0  =  x0. 
Entonces,  mediante  el  grado  de  convexidad  logarítmico,  vamos  a  comparar  las 
sucesiones  {xn}  e  {yn}. 


Teorema  1.11. 

Si  Lf(x)  >  Lg(x)  para  f(x)  >  0,  entonces  la  sucesión  {yn}  converge  a  s  más  rápi¬ 
damente  que  {xn}. 

Demostración  : 

En  nuestras  condiciones  {xn}  e  {yn}  son  sucesiones  decrecientes  a  s.  Por 
tanto,  mediante  inducción  vamos  a  probar  que  yn  <  xn  para  todo  ne  N  tal  que 
xn  *  s.  Además  yn  =  s  si  xn_-|AS  y  xn  =  s. 

Teniendo  en  cuenta  que  F(s)  =  G(s)  =  s  y  x0  *  s,  tenemos  que  x-|  -  y-j  = 

(  F  -  G  )(x0)  -  (  F  -  G  )(s)  y,  por  el  Teorema  del  Valor  Medio,  existe  e  (  s,  x0) 

tal  que  x-j  -  y-j  =  (F  -  G)’(^0)(x0-s).  Por  otra  parte,  como  (F  -  G)'(x)  =  Lf(x)  -  Lg(x) 
y  f(£0)  >  0  entonces  obtenemos  que  x-|  -  y-|  >  0. 

Supongamos  que  x^  >  y^  para  k  =  1 ,  2 . n-1 ,  si  xn  *  s  entonces  xn_i  a  s 

y  xn  -  yn  =  F(xn-l) '  G(yn-l)  -  (  F  ‘  G  )  (xn-1 )  puesto  que  G  es  una  función  cre¬ 

ciente  en  [s,b]  e  {yn}  c  [s,b].  Por  tanto,  procediendo  como  en  el  caso  k=1 ,  obte¬ 
nemos  que  se  verifica  xn  -  yn  >  0. 

Aplicando  los  argumentos  utilizados  anteriormente  es  sencillo  probar  que  si 
xn-1  *  s  y  xn  =  s  entonces  yn  =  s.  # 

Luego  hemos  probado,  mediante  e!  grado  de  convexidad  logarítmico,  que  a 
menor  convexidad  de  una  función  el  Método  de  Newton-Ftaphson  converge  más 
rápidamente  a  la  raíz  de  la  ecuación. 
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Ejemplo  1.4. 

Sean  las  funciones  f(x)  =  (  x3/216  )  -1 ,  g(x)  =  (  x2/36  )-1 . 

Veamos  que  ambas  tienen  las  mismas  características  en  [5,10],  además  de 
la  misma  raíz  s=6. 

f’(x)  =  x2/72  ;  f”(x)  =  x  /36  ;  Lf(x)  =  (2/3)  -  144/x3 

g’(x)  =  x/18  ;  g”(x)  =1/18  ;  Lg(x)  =  (1/2)  -  18/x2 

vemos  que  las  dos  son  convexas  y  crecientes  en  dicho  intervalo,  además  se  veri¬ 
fica  Lg(x)  <  Lf(x)  <=>  x3  +  108  x  -  864  >  0 


na» 

f 

g  1 

i 

10.000000000000000 

10.000000000000000  1 

2 

7.386666666666700 

¿ibbbbbbbbbbbbbb  ! 

3 

6.244023743014800 

6.047058823529500  ' 

4 

6.009412497423800 

6.000183108262900  ! 

5 

6.000014735026500 

6.000000002793900  ' 

6 

6.000000000036300 

¿.bbbbbbbbóbóbbbb  ¡ 

7 

6.000000000000100 

6.000000000000000  ' 

V _ / 


Tabla  5. 


y  por  tanto  Lg(x)  <  Lf(x)  en  [5,10],  luego  la  sucesión  {yn}  construida  con  la  función 
g  converge  a  la  raíz  más  rápidamente  que  la  {xn}  construida  con  f.  (  Tabla  5). 
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CAPÍTULO  II 

LA  CONVEXIDAD  Y  LA  CONCAVIDAD  EN  LA  CONVERGENCIA 
DE  PROCESOS  ITERATIVOS  PARA  LA  RESOLUCIÓN  DE 
ECUACIONES  NO  LINEALES 


2.0.  Introducción 

En  este  segundo  Capítulo  realizamos  un  análisis  de  la  convergencia  de  los 
dos  procesos  iterativos  de  tercer  orden  más  utilizados,  el  método  de  Chebyshev 
o  de  la  interpolación  cuadrática  inversa  [4],  [5],  y  el  método  de  Halley  [1],  [9]. 
Además,  a  partir  de  estos  procesos  iterativos  consideramos  una  familia  unipara- 
métrica  de  procesos  iterativos  de  tercer  orden. 

En  lo  que  sigue  consideraremos  f  e  c(P  ( [a,b] ),  p  >  4,  en  las  condiciones  de 
Fourier,  es  decir,  f(a)  f(b)  <  0,  f’(x)  *  0  y  f”  con  signo  constante  en  [a,b],  entonces 
existe  una  única  raíz  de 

(2.1.)  f(x)  =  0 

en  [a,b].  En  nuestro  estudio  podemos  suponer,  sin  pérdida  de  generalidad,  que  f 
es  una  función  convexa  y  estrictamente  creciente  en  [a,b],  en  otro  caso  basta 
cambiar  f(x)  por  f(-x),  -f(x)  o  -f(-x). 

A  lo  largo  de  todo  el  Capítulo  denotaremos  por  m(f)  =  mín.{f(x)l  x  e[a,b]}  y 
M(f)  =  máx  {f(x)l  x  e[a,b]}. 


2.1.  El  Método  de  Chebyshev. 

2.1.1.  Generalidades 

Dada  la  función  y  =  f(x),  denotamos  x  =  tp(y)  a  la  función  inversa  de  f.  Si  la 
aproximamos  en  un  punto  x0,  con  y0  =  f(x0),  entonces  cp(y)  ~  cp(yQ)  +  <p‘(yo)(y  -  yo) 
+  <P"(y0)(y  -  y0)2/2. 
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Hallar  la  raíz  de  la  ecuación  f(x)  =  0,  corresponde  a  encontrar  la  imagen  del 
cero  por  la  función  cp,  es  decir  9(0).  Entonces  <p(0)  =  <p(y0)  -  <p’(y0)  y0  + 
tp”(y0)(  y0)2/2,  por  tanto  tenemos  que 

(2.2.)  s  =  x0  -  <p’(y0)  f(x0)  +  cp”(y0)  f(x0)2/2. 

Hallemos  q>‘(y)  y  q>"(y),  para  poder  sustituir  su  valor  en  la  expresión  anterior. 

x  =  <p(y) 

y’=  f’(x),  y”  =  f ”(x) 


0'  (y)  =  —  =  -í—  =  — 
T  vy;  dy  dy  y' 
dx 


1 

f(x) 


r  d  dx-j/dXx  _  _d,dx_\  dx 
dx  dy  dy  _  dx^dy'dy 


jL(J_\í|x.  _  _  y_ 

dx  y'  dy  _  .3 


Sustituyendo  en  (2.2.),  tenemos 


s  ~  x 


f(X„)  l"(x0)f(x 


0  l'W  2f 


Entonces  consideramos 


'(x0) 

1  0  f'(Xo)  2  J 


Repitiendo  para  este  x-|  la  construcción  anterior  obtenemos  el  término  X2  y 
por  recurrencia  se  tiene 


f  (xn  i )  L  (xn1) 

(2.3.)  x  =  x  , - ^^[1  +  -i--"-1—  ] 

v  ’  n  n_  1  f  (x  )  2  J 

'  n- 1 ' 


que  es  el  proceso  iterativo  conocido  como  el  método  de  Chebyshev  [4],  [5], 
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2.1.2.  Discusión  de  la  convergencia 

A  continuación  vamos  a  estudiar  la  convergencia  de  la  sucesión  {xn},  dada 
por  el  proceso  iterativo  (2.3.). 


Teorema  2.1. 

Supongamos  que  Lf  >(x)  <  3  en  [a,b].  Si  f(x0)  >  0  entonces  la  sucesión  {xn}  es 
decreciente  y  converge  a  s. 

Además,  si  f(x0)  <  0  y  Lf  (x)  >  -2  en  [a,b]  entonces  la  sucesión  {xn}  es  cre¬ 
ciente  y  converge  a  s. 

Demostración: 

Consideramos  x0  e  [a,b]  tal  que  f(x0)  >  0  y  probaremos  que  xn>s  para  todo 
n  e  N. 

Como  f(x0)  >  0  se  sigue  que  x0  -  s  >  0  ,  por  hipótesis.  Además,  x-|  -s= 
F(x0)-F(s)=  F’(0O)  (x0-  s)  con  0O  e  (s,x0),  y  es  fácil  obtener  que  F’(x)  >  0  en  [a,b] 
si  Lf  >(x)  <  3,  teniendo  en  cuenta 


2  3-L  (x) 

(2.4.)  F'(x)  =  Lf  (x)  - 1 - 


Entonces  x-j  >  s  ,  y  por  inducción  xn  >  s  para  todo  n  e  N. 


Por  otra  parte,  xn-x  _  = 

n  n-1  f 


f(Xn-l)  Lf(Xn-l) 


'(XnJ 


[1  + 


]  <  0  pues 


f(xn)  >  0  para  todo  n  e  N.  Luego  {xn}  es  decreciente  y  pasando  al  límite  en  (2.3.), 
se  sigue  que  converge  a  la  única  raíz  s,  de  (2.1.). 

Si  f(x0)  <  0  la  demostración  es  análoga.  # 

Notar  que  cuando  consideramos  x0  tal  que  f(x0)  <  0,  la  condición  Lf  (x)  >  -2 
en  [a,b]  no  es  muy  restrictiva  puesto  que  Lf  (s)  =  0  y  Lf  (x)  >  0  en  [s,b]. 
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Definición  2.2. 

Sea  {xn}  una  sucesión  que  converge  a  s.  Diremos  que  {xn}  es  oscilante,  si 
x2  ( j  - 1  )  -  x2  j  -  s  >  x2j  + 1  >x2j.1  óbienx2j..,  >x2j  +  1  >s>x2j  >x2(j.1) 
para  j  >  0. 

En  los  teoremas  que  siguen  la  sucesión  {xn}  obtenida  va  a  ser  oscilante. 
Por  ello,  en  primer  lugar  vamos  a  obtener  una  condición  que  nos  asegure  que 

{xn>  c  [a,b]. 


Lema  2.3. 

Sea  8  =  máx.{-f(a),  3f(b)/2  },  y  supongamos  que  F  es  decreciente  y  Lf(x)  <  1 
en  [a,b].  Si  b  -  a  >  8/f'(a)  entonces  xne[a,b]  para  todo  neN. 

Demostración: 

Como  xn  =  F(xn_i)  y  F  es  una  función  decreciente  basta  probar  que  F(b)  >  a, 
F(a)  <  b  y  aplicar  inducción.  Al  ser  Lf(x)  <  1  en  [a,b]  entonces  F(b)  >  b  -  3f(b)/2f  '(b), 
y  a  partir  de  la  hipótesis  obtenemos  que  F(b)  >  a.  Análogamente  se  prueba  que 
F(a)  <  b,  teniendo  en  cuenta  que  Lf(a)  <0.  # 

Así,  podemos  suponer  que  el  intervalo  [a,b]  verifica  que  b-a  >  8/f'(a),  con  lo 
que  la  sucesión  {x^,-  [a,b].  Esta  condición  es  difícil  de  comprobar  en  la  práctica 
y  por  ello  en  cada  resultado  de  convergencia,  dando  condiciones  sobre  el  punto 
de  salida  aseguramos  que  {xn}  [a,b],  en  vez  de  exigir  que  el  intervalo  [a,b]  veri¬ 
fique  la  condición  dada  en  el  lema  anterior,  si  bien  en  determinadas  situaciones 
puede  ser  útil  dicha  condición. 


Teorema  2.4. 

Sea  xge[a,b]  con  f(xo)>0  y  xg  >  a  Si  Lf(x)e(3,5]  y  ILf(x)l<1  en  [a,b], 

entonces  la  sucesión  {xn}  es  oscilante  y  converge  a  s. 

Demostración: 

En  primer  lugar,  veamos  que  mediante  las  condiciones  dadas  para  el  punto 
inicial  xg,  un  procedimiento  inductivo  prueba  que  {xn}  está  contenida  en  [a,b].  Es 
claro  por  hipótesis  que  x-|=F(xg)<xo<b.  Por  otra  parte, 

Xl  =  F(x0)  >  x0  -  3f(x0)/2f(x0)  >  x0  -  3f(b)/2f’(b)  >  a 
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por  ser  0<Lf(x0)<1  y  f/f’  creciente  en  [a,b],  Por  tanto,  x-|e[a,b]. 

En  las  hipótesis  del  teorema  y  teniendo  en  cuenta  (2.4.),  se  sigue  fácilmente 
que  F’(x)<0  y  IF’(x)l<1  en  [a,b].  Si  f(xg)>0,  es  decir  xg>s,  tenemos  que 

xrs  =  F(xg)-F(s)  =  F’(^g)-(x0-s)  con  ^ge(s,x0), 

luego  x-|<s.  Supongamos  que  X2k_i<s  y  X2k>s  para  k=1,2,...,n-1 .  De  manera  aná¬ 
loga  al  caso  k=1 ,  como  F’(x)<0,  se  prueba  que  x2n.i<s  y  X2n^s.  Por  lo  tanto,  {xn} 
es  una  sucesión  oscilante. 

Además,  como  IF’(x)l<1  para  todo  xe[a,b],  se  sigue  que 

lx-| -si  <  M(IF’(x)l)  IXg-Sl  <  Ixg-sl. 

Supongamos  a  continuación  que 

Ix^-sl  <  M(IF’(x)l)k|xQ-sl  para  k=1,2,...,n-1. 

De  manera  análoga  al  caso  k=1 ,  tenemos 

lxn-sl  <  M(IF’(x)l)-lxn.-|-sl  <  M(IF’(x)l)nlxg-sl. 

Como  M(IF’(x)l)<1 ,  entonces  {xn}c  [a,b]  y  lím  lxn-sl=0.  Por  lo  tanto,  {xn}  con¬ 
verge  as.  # 

Teorema  2.5. 

Sea  xge  [a,b]  con  f(x0)>0  y  x0  >  a  +  Si  Lf(x)>5  y  ILf(x)l  < 

[a,b],  entonces  la  sucesión  {xn}  es  oscilante  y  converge  a  s. 

Demostración: 

En  primer  lugar,  se  prueba  como  en  el  teorema  anterior  que  x-|e[a,b]. 

Por  otra  parte,  es  inmediato  comprobar  que  IF’(x)l<1  en  [a,b],  teniendo  en 
o  2 

cuenta  (2.4.)  y  que  Lf(x)¿  <  en  [a,b].  El  resto  de  la  demostración  se 

sigue  de  forma  análoga  a  la  del  teorema  anterior.  # 

fía) 

En  el  caso  en  que  xge  [a,b]  con  f(xg)<0,  basta  con  ver  que  xg  >  b  +  para 
el  mismo  resultado  que  en  los  teoremas  anteriores.  Notemos  que  en  ambas  situa¬ 
ciones  tenemos  que  asegurar  la  existencia  de  un  punto  xge[a,b]  verificando  la 

condición  xg  >  a  es  decir,  en  los  teoremas  anteriores  ha  de  verificarse  que 

3fíbt 

a  +|f(b)  -  b-  Ahora  bien,  como  en  la  práctica,  habitualmente,  el  dominio  de  f  es 
todo  R,  esta  cantidad  nos  permitirá  fijar  el  intervalo  [a,b]. 
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c 


2. 1.3.  Notas  prácticas 

A  continuación,  pasamos  a  estudiar  qué  ocurre  en  la  práctica.  Es  decir,  dada 
la  ecuación  (2.1),  ¿siempre  podemos  aproximar  su  raíz  mediante  el  método  de 
Chebyshev?.  Veamos  que  en  la  práctica  esta  pregunta  tiene  respuesta  afirmativa. 
Según  los  valores  de  m(Lf)  y  M(Lf>),  se  pueden  presentar  las  siguientes  situacio¬ 
nes 


i)  M(Lf-)<3 

ii)  M(Lf>)^5  y  m(Lp)>3 

iii)  3<M(Lf’)<5  y  m(Lf’)<3 

iv)  M(Lp)>5  y  m(Lf’)<5 

v)  m(Lf’)>5, 


como  puede  verse  en  la  Fig.27. 


Notar  que  en  las  situaciones  (i),  (ii)  y  (v)  la  convergencia  de  la  sucesión  {xn} 
está  asegurada  por  los  Teoremas  2.1.,  2.4.  y  2.5.  respectivamente,  con  algunas 
condiciones  sobre  Lf . 

Por  otra  parte,  en  los  casos  (iii)  y  (iv)  no  sabemos  nada  acerca  del  creci¬ 
miento  del  operador  F.  En  estos  casos,  sólo  podemos  asegurar  que  la  sucesión 


92 


LA  CONVEX.V  LA  CONCA.  EN  LA  CONVERGENCIA  DE  PROCESOS  ITERATIVOS  PARA  LA  RESOLUCIÓN  DE  ECUACIONES  NO  LINEALES 


{xn}  converge  a  s.  Así  al  igual  que  en  el  caso  en  que  {xn}  es  una  sucesión  osci¬ 
lante,  debemos  asegurar  que  {xn}  está  contenida  en  [a,b],  Aquí,  se  pueden  hacer 
dos  observaciones.  La  primera  de  ellas  es  la  misma  que  la  que  sigue  al  Lema  2.3.; 
y  la  segunda,  consiste  en  encontrar  un  intervalo  en  el  cual  cualquier  punto  perte¬ 
neciente  a  él,  pueda  ser  punto  inicial  para  la  sucesión  {xn}.  El  Lema  siguiente 
hace  referencia  a  esta  última  consideración. 

Lema  2.6. 

Supongamos  que  IF’(x)l<1  y  ILf(x)l<1  en  [a,b], 

i)  Sea  ote  [a,b]  tal  que  f(a)<0  con  <  a-a.  Si 

x0e[a,b]  con  f(x0)  >  0  entonces  {xn}cz  [a,b]. 

ii)  Sea  pe  [a,b]  tal  que  f(P)>0  con  >  p-b.  Si  x0e  [a,b]  con  f(xg)  <  0  enton¬ 
ces  {xn}c  [a,b]. 

Demostración: 

Para  probar  (i),  veamos  en  principio  que  x-|e[a,b].  Es  claro  que  x-j<xg  y 

v  -  FÍX  W  v  3f(x0>  >V  3í(b)  >  a  _ 

X1  -  F(xo)  >  x0  '2f’(x0)  “  Xn  "  2f’(b)  ~  x0'a  +  a  >  a- 
por  ser  Lf(x)<1  y  f/f’  creciente  en  [a,b].  Luego,  x-|e[a,xo)  c  [a,b].  A  continuación, 
distinguimos  dos  casos:  x-|<s  y  x-|>s.  En  primer  lugar,  si  x-|<s,  entonces 
-1  <Lf(x-|  )<0  y  como  f(x-|  )<0  se  sigue  que  x-|<x2.  Además,  se  tiene  que 
lx2-sl<lxi -sl<lx0-sl  al  ser  IF’(x)l<1.  Por  tanto,  x2e(x-|,xo)c  [a,b].  Ahora,  aplican¬ 
do  hipótesis  de  inducción  sobre  n,  obtenemos  que  xne(xi,xg)c:  [a,b]  para  todo 
neN.  En  segundo  lugar,  si  x-ps,  entonces  s<xi<xg.  Por  tanto,  x-|e[a,b]  y 

x2  =  F(x-, )  >  X-,  -  >  x-,  -  J¡^  >  xra+a  >  a. 

Luego,  x2e[a,x-|)c  [a,x0)c  [a,b]. 

Ahora  distinguimos  de  nuevo  dos  casos.  Por  una  parte,  si  existe  algún 

xj<e(s,b]  con  k=1 n,  se  prueba  análogamente  al  caso  anterior  que 

xn+-|e[a,xn)  <z  [a,x0)  c  [a,b],  luego  {xn}  c  [a,b].  Por  otra  parte,  si  existe  algún 
xn  e  [a,s)  con  x^  e  (s,b]  para  k  =  1,...,n-1,  se  prueba  de  la  misma  forma  que  el 
caso  anterior  que: 

xn+1e(xn-xn-l)c  (xn-xo)c  ta-bl>  luego  {xn}  c  [a,b]. 

El  caso  (ii)  se  sigue  de  forma  análoga.  # 

En  el  teorema  siguiente,  elegimos  xg,  a  y  p  en  [a,b]  verificando  las  condi¬ 
ciones  del  Lema  2.6. 
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Teorema  2.7. 

Supongamos  que  M(ILfl)  <  '\[^y  K=máx.{l3-M(Lf’)l,l3-m(Lf')l)}. 

i)  Si  3<M(Lf)<5,  m(Lf’)<3  y  ILf(x)l<1  en  [a,b],  entonces  la  sucesión  {xn} 

converge  a  s. 

ii)  Si  M(Lf)^5,  y  m(Lf’)<5,  entonces  la  sucesión  {xn}  converge  a  s. 

Demostración: 

En  el  caso  (i),  veamos  primero  que  {xn}  está  contenida  en  [a,b].  Para  ello, 
basta  probar  que  IF’(x)l<1  en  [a,b]  y  aplicar  el  Lema  2.6.  En  efecto, 
IF’(x)l  <  2  ILf(x)2|<1  en  [a,b]  por  hipótesis.  Por  otra  parte, 

Ix-j -si  <  M(IF’(x)l)  lx0-sl  <  Ixq-sI. 

Supongamos  ahora  que  IXp-sl  <  M(IF’(x)l)P|x0-sl  para  p=1,2 . n-1. 

Así,  de  forma  análoga  al  caso  p=1 ,  obtenemos 

lxn-sl  <  M(IF(x)l)lxn_rsl  <  M(IF’(x)l)nlx0-sl  <  Ixq-sI. 

Por  lo  tanto,  como  M(IF’(x)l)<1 ,  tenemos  que  lím  lxn-sl=0.  Entonces,  {xn}  con¬ 
verge  a  s. 

Para  probar  (ii),  basta  tener  en  cuenta  que,  en  las  hipótesis  anteriores,  K>2  y 
por  tanto  ILf(x)l<1  en  [a,b].  La  convergencia  de  {xn}  se  prueba  como  en  el  caso 
anterior.  # 

A  partir  de  este  resultado,  podemos  deducir  que  siempre  podemos  resolver 
la  ecuación  (2.1),  aplicando  el  método  de  Chebyshev  para  todos  los  valores  de 
M(Lp)  y  m(Lf>)  con  algunas  restricciones  sobre  los  valores  de  Lf(x)  en  [a,b],  Pero 
teniendo  en  cuenta  que  Lf(s)=0,  es  claro  que  en  un  entorno  de  la  raíz  s,  siempre 
podemos  aplicar  estos  resultados. 


Ejemplo  2.1. 

Dada  la  función  f(x)=x3-x2-x-2,  veamos  que  el  método  de  Chebyshev  es  con¬ 
vergente  a  la  raíz  de  la  ecuación  f(x)=0  .  Es  fácil  obtener  que 

6(3x2-2x-1) 


L,(x)  =  a3:^M6x-2_)y 


Lf(x)  = 


(3x2-2x-1)2  7  “IX'V  (6x-2)2 

Entonces,  se  sigue  de  forma  sencilla  que  Lf(x)<1/2  en  R+  y  2/3>Lf(x)>-2  para 

todo  x  >1’51763.  Por  tanto,  si  restringimos  el  intervalo  inicial  al  intervalo  [1’52,4], 
y  tomamos  por  ejemplo  xq=1’6,  entonces  por  el  Teorema  2.2.  la  sucesión  {xn} 
dada  por  (2.3.)  es  creciente  y  converge  a  la  solución  de  f(x)=0,  ver  Tabla  6. 


iter. 

xn 

0 

1 .600000000000000 

i 

1 .808984924897836 

2 

1.989904870759261 

3 

1 .999999075443373 

4 

2.000000000000000 

Tabla  6 


94 


LA  CONVEX.V  LA  CONCA.  EN  LA  CONVERGENCIA  DE  PROCESOS  ITERATIVOS  PARA  LA  RESOLUCIÓN  DE  ECUACIONES  NO  LINEALES 


2.2.  El  método  de  Halley. 

2.2.1.  Generalidades 

Es  conocido  [12]  que  dada  una  ecuación  no  lineal  (2.1.),  si  para  su  resolución 
aplicamos  el  método  de  Newton,  también  llamado  método  de  la  tangente,  este 
proceso  iterativo  posee  convergencia  cuadrática  con  constante  del  error  asintóti- 
co  lf”(s)/2!  I,  siendo  s  la  raíz  de  dicha  ecuación.  Entonces,  un  procedimiento  para 
aumentar  la  velocidad  de  convergencia  de  este  proceso  iterativo,  consiste  en  con¬ 
seguir  que  esta  constante  del  error  asintótico  sea  cero,  con  lo  cual  el  proceso  ite¬ 
rativo  pasaría  a  tener,  al  menos,  convergencia  cúbica.  Para  ello,  en  general,  dada 
una  función  arbitraria  g  sabemos  que  el  producto  f(x)  g(x)  preserva  los  ceros  de 
la  función  dada  f.  Por  ello  vamos  a  resolver  la  nueva  ecuación  no  lineal  f(x)g(x)  = 
0.  Así,  podemos  plantearnos  cómo  ha  de  ser  la  función  g  para  que  se  verifique 
(fg)”(s)  =  0,  con  lo  cual  ya  obtendríamos,  al  aplicar  el  método  de  Newton,  un  pro¬ 
ceso  iterativo  con  convergencia  cúbica.  Para  ello,  en  un  entorno  apropiado  [a,b] 
de  la  raíz  s,  es  suficiente  suponer  que  f,ge  c(  2([a,b])  y  que  f’(x)g(x)  *  0.  Notemos 
que  f(s)  =  0  y  además  (fg)”(s)  =  0.  Entonces,  desarrollando  esta  segunda  deriva¬ 
da  podemos  escribir 


(2.5.) 


g'(s) 

g(s) 


f"(s) 

2  f'(s) 


Así,  integrando  (2.5.)  sobre  el  intervalo  [a,b]  eligiendo  apropiadamente  la 
constante  de  integración,  obtenemos  que  la  función  g  =  (f’)‘ 1/2  verifica  la  condi¬ 
ción  (2.5.).  Entonces  aplicando  el  método  de  Newton  a  la  ecuación  f(x)g(x)  =  0 
obtenemos  el  siguiente  proceso  iterativo 


*n  = 


n- 1 


[f,(xn.1)/f(xn.1)]-  [f"(xn.1)/2f,(xn.1)] 


para  n  >  1 . 


Mediante  el  grado  de  convexidad  logarítmico  este  proceso  iterativo,  denomi¬ 
nado  método  de  Halley,  [9]  admite  la  expresión 


(2.6.)  x  =  F(x  ) 

n  n- 1 


con 


F(x) 


f  (x)  2 

f  (x)  2  -  L  (x) 
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2.2.2.  Discusión  de  la  Convergencia 

Dedicamos  esta  Sección  a  estudiar  la  convergencia  del  método  de  Halley  a 
partir  del  grado  de  convexidad  logarítmico  de  las  funciones  f  y  f.  La  expresión  de 
este  método  viene  dada  por  (2.6).  A  continuación,  estudiamos  la  convergencia 
para  el  método  de  Halley  en  el  caso  en  que  la  sucesión  {xn}  sea  monótona  a  la 
raíz. 


Teorema  2.8. 

Supongamos  Lf  >  (x)<3/2  en  [a,b].  Si  f(xg)>0  y  Lf(x)<2  en  [a,b],  entonces  la 
sucesión  {xn}  dada  por  (2.6)  es  decreciente  y  convergente  a  s.  Si  f(xg)<0,  enton¬ 
ces  {xn}  es  creciente  y  convergente  a  s. 

Demostración: 

Supongamos  que  xg  e  [a,b]  tal  que  f(xg  )>0.  Veamos  que  xn>s  y  xn .  i  >xn  para 
todo  neN.  En  primer  lugar,  xg-s>0  por  ser  f(xg)>0.  Además, 

xi-s  =  F(xg)-F(s)  =  F'fég)-(xg-s)  con  Sg[(s,x0) 

Ahora  bien,  teniendo  en  cuenta 


Lf(x)2 

(2.7)  F'(x)  = - - -  [3  -  2L  (x)  ] 

[2  -  Lf(x)] 

es  claro  que  F’(x)>0  en  [a,b]  por  ser  Lf  -  (x)<3/2.  Luego  x-j  >s.  Por  otra  parte,  como 

Lf  (x)<2,  tenemos  que  x-|  -xg  =  - ? - <  0.  A  continuación,  aplicando 

f’  (x0)  2-Lf(x0) 

hipótesis  de  inducción  sobre  n,  se  sigue  de  forma  fácil  que  xn>s  y  xn .  i>xn  para 
todo  ne  N.  Entonces,  la  sucesión  {xn}  es  decreciente  y  acotada.  Por  lo  tanto,  {xn} 
es  convergente.  Pasando  al  límite  en  (2.6.)  obtenemos  que  precisamente  este 
límite  ha  de  ser  la  única  raíz  s  de  f(x)  =  0.  # 

En  los  teoremas  siguientes  la  sucesión  {xn}  dada  por  (2.6.)  no  va  ser  monó¬ 
tona.  Por  eso  vamos  a  obtener  en  primer  lugar  condiciones  que  nos  aseguren  que 
{xn}c[a,b]. 


Lema  2.9 

Supongamos  I  F'(x)  I  <  1  y  I  Lf(x)l  <  1  en  [a,b]. 
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(i)  Sean  x0,  a  e  (a,b)  tales  que  f(oc)  <  0  y  f(x0)  >  O,  si  2f(x0)/f'(x0)  <  a  -  a 
entonces  {xn}  c  [a,b]. 

(ii)  Sean  x0,  p  e  (a,b)  tales  que  f(P)  >  O  y  f(x0)  <  O,  si  f(x0)/f'(x0)  >  -b  +  p 
entonces  {xn}  c  [a,b]. 

Demostración: 

(i)  Como  x-|  =  F(x0)  >  x0  -  2  f(x0)/f'(x0)  >  x0  -  a  +  a  >  a,  es  claro  que  x-|  e 
(a,  x0)  c=  [a,b] .  Distinguimos  dos  casos  x-|  <s  y  x-|  >s. 

(a)  Si  x-j  <  s  entonces  -1  <  Lf  (x-j )  <  O  y  como  f(x-| )  <  O  se  sigue  que  x-|  < 
X2  Además,  se  tiene  que  I  X2  -  s  I  <  I  x-|  -  s  l<  I  xQ  -si  al  ser  IF’(x)l  <  1.  Se 
deduce  por  tanto  que  X2  e  (x-| ,  x0)  c  [a,b] .  Un  proceso  inductivo  prueba  que 
xn  e  (x-| ,  x0)  c  [a,b], 

(b)  Cuando  x-|  >  s  entonces  s  <  x^  <  xQ,  por  tanto  x-|  e[a,b]  y  *2  =  F(x1 )  - 
x-|  -  2  f(x1  )/f’(x-| )  >  x-|  -  2  f(x0)/f’(x0)  >  x0  -  a  +  a  >  a  puesto  que  f/f'  es  una  fun¬ 
ción  creciente  en  (s,  b).  Luego,  x2e(a,Xf )  c  (a,  x0)  c  [a,b] .  De  nuevo  podemos 
distinguir  dos  casos: 

(b.1)  Si  xk  e  (s,b)  para  todo  k  =1,.,n-1,  podemos  probar  de  manera  análoga 
que  xn  e  (  a,  x0  )  y  por  tanto  {xn}  <=  [a,b]. 

(b.2)  Si  existe  algún  xk  <  s  con  x-| . xk-  1e  (s,  b)  se  puede  probar,  como 

en  la  parte  (a)  que  xk+i  e  (  xk,  x0  )  c  [a,b] .  Por  tanto  {xn}  c  [a,b]. 

(ii)  Se  prueba  de  la  misma  forma.  # 

Las  condiciones  del  resultado  anterior  pueden  suavizarse  en  el  caso  de  que 
podamos  afirmar  que  {xn}  es  oscilante,  así  tenemos  que  para  Lf  (x)>3/2  en  [a,b], 
la  sucesión  {xn}  dada  por  (2.6.)  es  oscilante,  y  por  tanto,  tenemos  que  asegurar 
que  x-|e[a,b],  lo  que  nos  lleva  a  considerar  únicamente  una  condición  sobre  el 
punto  de  salida. 


Teorema  2.10. 

Sea  x0e  [a,b]  con  f(xo)>0  y  x0  >  a  +  Si  Lf  >(x)e  (3/2,2]  y  ILf  (x)l<1  en  [a,b], 
entonces  la  sucesión  {xn}  converge  a  s.  Además,  esta  sucesión  es  oscilante. 
Demostración: 

Teniendo  en  cuenta  la  condición  dada  para  el  punto  de  salida  x0  y  que  f/f’  es 
creciente  en  [a,b],  se  deduce  fácilmente  que  x-|  e[a,b]. 
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A  continuación,  veamos  que  la  sucesión  {xn}  es  oscilante  y  está  contenida 
en  [a,b],  Si  f(xo)>0,  es  decir,  x0>s,  entonces  x-|-s  =  F(x0)-F(s)=  F’(^0)(x0-s)  con 
¡;0e(s,x0),  luego  x-(  <  s  y  lx-|-sl  <  M(IF’l)lx0-sl.  Ahora,  es  fácil  probar  por  induc¬ 
ción  que  x2k  >  s,  X2|<+i<s  y  lx|<+-|-sl  <  M(IF’l)k+1  lx0-sl  para  todo  k>0.  Luego, 
{xn}  es  oscilante.  Además,  como  M(IF’I)<1  en  [a,b],  obtenemos  que  la  sucesión 
{xn}  c  [a,b],  la  sucesión  {lxn-sl}  es  decreciente  y  converge  a  cero.  Por  tanto  {xn} 
converge  as.  # 


Teorema  2.11. 

Sea  x0e  [a,b]  con  f(xo)>0  y  x0  >  a  +  Si  Lf-(x)>2  y  ILf(x)l  <  ^(L¡!)f-^  en 

[a,b],  entonces  la  sucesión  {xn}  converge  a  s  y  es  oscilante. 

Demostración: 

De  la  misma  forma  que  en  el  Teorema  anterior  se  sigue  que  x-|  e  [a,b].  Por  otra 
parte,  se  comprueba  fácilmente  que  IF’(x)l<1  en  [a,b].  Basta  tener  en  cuenta  (2.7.) 
y  que  2(M(Lf>)-2)Lf(x)2  +  4Lf(x)-4<0  por  ser  Lf(x)  <  ~1+  en  [a,b].  El  resto 

de  la  demostración  es  similar  a  la  del  teorema  anterior.  # 

En  el  caso  en  que  x0e  [a,b]  con  f(xo)<0  y  x0  >  b  +  f(a)/f’(a),  se  siguen  los  mis¬ 
mos  resultados  en  los  dos  Teoremas  anteriores.  Notemos  que  en  ambas  situa¬ 
ciones  tenemos  que  asegurar  la  existencia  de  un  punto  x0e[a,b]  verificando  la 
condición  correspondiente,  es  decir,  por  ejemplo  en  el  caso  de  los  teoremas  ante¬ 
riores  ha  de  verificarse  que  a  +  <  b.  Ahora  bien,  como  en  la  práctica,  habi¬ 

tualmente,  el  dominio  de  f  es  todo  R,  esta  cantidad  nos  permitirá  fijar  el  intervalo  [a,b]. 

2.2.3.  Notas  prácticas 

Probamos  a  continuación  que  en  la  práctica  siempre  se  puede  aplicar  el 
método  de  Halley  para  resolver  la  ecuación  (2.1.).  Se  obtienen,  a  partir  de  los 
valores  de  M(Lf)  y  m(Lf>),  los  siguientes  casos 

i)  M(Lf  )<3/2 

ii)  M(Lf)<2  y  m(Lf>)>3/2 

iii)  3/2<M(Lf-)<2  y  m(Lf>)<3/2 

iv)  M(Lf)>2  y  m(Lf’)^2 

v)  m(Lf>)>2, 
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como  puede  verse  en  la  Fig.  28. 


Notar  que  en  los  casos  (i),  (i¡)  y  (v)  la  convergencia  de  {xn}  está  estudiada  en 
los  Teoremas  2.8.,  2.10.  y  2.1 1 .  respectivamente  con  algunas  condiciones  sobre  Lf. 

A  continuación,  estudiamos  los  casos  (iii)  y  (iv).  En  ellos  no  sabemos 
cómo  va  a  ser  la  convergencia  de  la  sucesión  {xn},  por  ello  para  asegurar  que 
{xn}c  [a, b]  podemos  aplicar  el  Lema  2.9.,  en  el  que  se  asegura  este  hecho  dando 
condiciones  sobre  el  punto  de  salida  o  bien  aplicar  el  siguiente  resultado  donde 
sólo  se  fija  la  ubicación  del  punto  de  salida  respecto  de  la  raíz  de  (2.1 .). 


Lema  2.12. 

Supongamos  que  IF’(x)l<1  y  ILf(x)l<1  en  [a,b]. 

i)  Sean  xg,ae[a,b]  tales  que  f(a)<0  y  f(xg)>0.  Si  ^rff$  <  a-a,  entonces 
{xn}  c  [a,b], 

ii)  Sean  XQ,Pe[a,b]  tales  que  f(p)>0  y  f(xg)<0.  Si  >  p-b,  entonces' 
{xn}  c  [a,b]. 
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Demostración: 

Totalmente  análoga  a  la  del  Lema  2.9.  # 

En  el  teorema  siguiente,  elegimos  xq,  oc  y  p  en  [a,b]  verificando  las  condicio¬ 
nes  del  Lema  2.12. 


Teorema  2.13. 


Sea  K  =  máx.{l3-2M(Lf)l,  l3-2m(Lf-)l}.  Supongamos  que  M(ILfl)  <  —  ^  . 
Entonces, 

i)  Si  3/2<M(Lf>)<  2,  m(Lf’)<3/2  y  ILf(x)l<1  en  [a,b],  la  sucesión  {xn}  converge  a  s. 
¡i)  Si  M(Lf>)>2,  y  m(Lf)^2,  la  sucesión  {xn}  converge  a  s. 

Demostración: 

Para  el  caso  (i),  teniendo  en  cuenta  (2.7.)  y  que  K  =  máx.{l3-2M(Lf)l,  13- 
2m(Lf’)l}  en  [a,b],  se  sigue  que 

lh  W  <(2-M(ILfl))2' 


p 

Luego,  IF’(x)l<1  en  [a,b]  por  ser  M(ILfl)  <  Ahora,  aplicando  el 

Lema  2.12.,  es  fácil  comprobar  que  la  sucesión  {xn}  c  [a,b]. 

En  el  caso  (ii),  basta  tener  en  cuenta  en  las  hipótesis  del  teorema  que  K>1 , 
y  por  tanto  M(ILfl)<1 .  Luego,  ILf(x)l<1  en  [a,b].  Por  otra  parte,  debemos  probar  que 
IF’(x)l<1  en  [a,bj,  es  decir,  (K-l)M(ILfl)2  +  4M(ILfl)  -  4  <  0.  Condicción  que  se  sigue 
fácilmente  por  ser  0  <  M(ILfl)  <1+ 

El  resto  de  la  demostración,  para  ambos  casos,  es  análoga  a  la  del  Teorema 
2.10.,  si  bien  no  podemos  asegurar  que  {xn}  sea  oscilante.  # 

Se  sigue  de  este  estudio  que  dada  la  ecuación  (2.1.)  siempre  podemos  apli¬ 
car  el  método  de  Halley  para  resolverla,  teniendo  en  cuenta  algunas  restricciones 
para  Lf(x),  que  en  realidad,  al  ser  Lf(s)=0,  siempre  podemos  encontrar  un  entor¬ 
no  de  s  en  el  cual  se  pueden  aplicar  estos  resultados. 


Ejemplo  2.2. 

Veamos  que  el  método  de  Halley  es  convergente  a  la  raíz  de  la  ecuación 
f(x)=-8+ex  '4=0  en  el  intervalo  [1,5].  Se  sigue  de  forma  sencilla  que 
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Lf(x)  - 


(1+2x2)(ex2~4) 

2x2ex2"4 


y  Lf-(x)  = 


2x2(3+2x2) 

(1+2x2)2 


Entonces,  es  fácil  probar  que  Lf>(x)<3/2  en  r+  y  Lf(x)<2  para  todo  x>0.  Por 
tanto,  podemos  tomar  como  punto  inicial  xq  para  (2.6)  cualquier  punto  del  inter¬ 
valo  [1 ,5]  que  verifique  f(xo)>0.  Por  ejemplo,  xg=4.  Así,  obtenemos  por  el  Teorema 
2.6.  que  la  sucesión  {xn}  dada  por  (2.6.)  es  decreciente  y  convergente  a  la  solu¬ 
ción  de  f(x)=0,  ver  Tabla  7. 


iter. 

xn 

0 

4.000000000000000 

1 

3.741961670527600 

2 

3.465034280857143 

3 

3.165567038560662 

4 

2.846385484142443 

5 

2.563813967367478 

6 

2.467671172927391 

7 

2.465652373241162 

8 

2.465652356209171 

9 

2.465652356209171 

Tabla  7 


2.2.4.  Apéndice 

Como  ya  se  ha  indicado  en  la  construcción  del  método  de  Halley,  éste  se 
puede  obtener  aplicando  el  método  de  Newton  a  la  función 


En  lo  que  sigue,  supondremos  f  en  las  condiciones  habituales.  Estudiando  la 
convergencia  del  método  de  Newton  para  la  función  h,  mediante  el  grado  de  con¬ 
vexidad  logarítmico,  obtenemos  un  resultado  de  convergencia  global  para  el 
método  de  Halley,  utilizando  para  ello  únicamente  resultados  de  convergencia 
para  el  método  de  Newton.  Si  bien  son  resultados  ya  obtenidos  anteriormente  tie¬ 
nen  interés  en  cuanto  a  la  técnica  utilizada. 
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Teorema  2.1 4V 

(i)  Si  Lf’(x)  <  3/2  en  [a,b]  entonces  {xn}  es  una  sucesión  decreciente  que  con¬ 
verge  a  s. 

(ii)  Si  Lf(x)  e  (  3/2,  2  )  y  Lf(x)  <  1  en  [a,b],  para  x0  >  a  +  2  f(b)/f'(a),  enton¬ 
ces  la  sucesión  {xn}  es  oscilante  y  converge  a  s. 

Demostración  : 

(i)  Es  inmediato  que  h  tiene  un  punto  de  inflexión  en  s,  ya  que 


h"(x)=  - f'(x)  ;  L  (x)  [  3  -  2  L  (x)  ]  . 

4  [  f  (x)  ] 

Si  Lf(x)  <  3/2  ,  se  sigue  que  h  es  una  función  cóncava  en  (a,s)  y  convexa  en 
(s,b)  .  Por  otra  parte  ,  como  h”  es  una  función  positiva  en  (s,b),  la  función 
h’(x)=f’(x)'l/2(2-Lf(x))/2  es  creciente  en  (s,b).  Además,  h’(s)  >  0,  y  por  tanto  h  es 
una  función  creciente  en  (s,b).  Aplicando  el  método  de  Newton  a  la  función  h, 
obtenemos  una  sucesión  decreciente  {xn}  cuyo  límite  es  s. 

(ii)  Como  F(x)  =  x  -  h(x)/h’(x) ,  se  sigue  de  (2.6.)  que 


Lh(x)  = 


_d_  r  x  h(x)  , 

dx  1  h'(x)  J 


f(x) 

f  (X) 


C(x)  ] 


con 


C(x)  = 


2 

2  -  L(x)  ' 


Como  es  una  función  negativa  en  [a,b],  obtenemos  que 


LJx)  >-1  «  2  -  (1-L(x)) 


2 _  f  (x)  2  Lf(x) 


2  -  L(x)  f  '(x)  [  2  -  L  (x)  f 


>0  . 


Entonces,  I  L^x)  I  <  1  <=>  f’(x)  (  2  -  Lf(x)  )  -  f(x)  L!f(x)  >  0. 


Por  otra  parte  ,  L'  (x)  =  [  1  +  L  (x)  (  L  (x)  -  2  )  ]  . 

'  f  (x)  '  ’ 


Por  tanto,  obtenemos  que  I  L^(x)  I  <  1  sí  y  sólo  sí  2-Lf(x)  > 
Lf(x)  [  1  +  Lf(x)  (Lp(x)  -  2 )  ].  Teniendo  en  cuenta  que  Lp(x)e  (3/2,2)  y  Lf(x)  <  1  en  [a,b], 
se  sigue  que  IL^(x)l  <1  en  [a,b].  Así,  existe  M  €  (0,1)  tal  que  I  L^(x)  I  <  M  en  [a,b]. 
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Notemos  que  F(x0)  e  (a,s).  Por  otra  parte,  como  x-|  -  s  =  F(x0)  -  F(s)  = 
Ffco)  (xo  '  s)  Para  ^o  6  (s’xo)  V  F(x)  =  LhM-  se  si9ue  que  I  x-|  -  s  I  <  M  I  x0  -  s  I. 
Por  inducción  ,  obtenemos  que  F(x2n)  e  (s,b)  y  F(x2n+-|)  e  (a,s)  para  n>0  ,  con 
lo  que  la  sucesión  es  oscilante. 

Además  ,  xn  -  s  =  F(xn_i)  -  F(s)  =  F'(£n--|)  (xn_i  -  s)  para  £n_i  un  punto  del 
intervalo  (mín.fx^  s},  máx.jx^  s}  ).  Entonces  lxn  -  si  <  M  lxn_-j  -  si.  Reiterando 
este  proceso,  obtenemos  que  lxn  -  s  I  <  Mn  I  xQ  -si  y  por  tanto  límn  xn  =  s.  # 


2.3.  Familia  de  procesos  iterativos  de  tercer  orden 


2.3. 1.  Discusión  de  la  Convergencia 

A  continuación,  vamos  a  considerar  una  familia  de  procesos  iterativos  del  tipo 
Chebyshev-Halley  [1 0]  con  al  menos  convergencia  cúbica  y,  mediante  los  valores 
de  Lf>,  obtendremos  un  resultado  de  convergencia  global  para  todos  ellos.  Para 
terminar  este  párrafo,  por  medio  de  un  estudio  de  optimización,  construiremos  un 
procedimiento  práctico  para  elegir  el  proceso  más  rápido  de  la  familia  dada  para 
resolver  (2.1 .). 

Por  otra  parte,  es  conocido  [5]  que  un  proceso  iterativo  dado  por  la  expresión 


f(x  J 

'  n-r 


xn=x  .  .  ..  ■  H(L.(X  )) 

n  n-1  f  (x  J  f  n-1" 

'  n-r 


donde  la  función  H  verifica  que  H(0)  =1,  H"(0)  =  1/2  y  IH"(0)  I  <  +  «,  tiene  al 
menos  convergencia  cúbica.  Consideramos  la  familia  de  procesos  iterativos 


(2.8) 


f  (x  ) 

xn,a  =  Fa(xn.1ia)=xn.1ia-  H(a’Lf(Xn-1,a) 

'  n- 1 ,  o! 


con 


H(a,x)  =1  +  ——  para  a  e  R  *  =  R  -  {0}  . 

2-  — 
a 

Es  fácil  ver  que  todos  ellos  tienen  al  menos  convergencia  cúbica.  Notar  ade¬ 
más  que  para  a  =  1  obtenemos  (2.6.)  y  si  consideramos  F±  ^  =  I  i  m  Fm 

.  m  — >  ±  oo 

obtenemos  (2.3.). 
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A  continuación  damos  un  resultado  de  convergencia  global  para  los  procesos 
iterativos  de  esta  familia. 


Teorema  2.15. 


Sea  a  eR*  y  x0>a  e  [a,b]  con  f(x0  a  )  >  0. 

i)  Si  a  <  0  y  Lf-(x)  <  3  -  3/(2a)  en  [a,b],  entonces  la  sucesión  {x0  a}  es  decre¬ 
ciente  a  s. 

ii)  Si  0  <  a  <  1 ,  Lf(x)  <  2  a  y  Lf -(x)  <  3  -  3/(2a)  en  [a,b],  entonces  la  sucesión 
{xn  a}  es  decreciente  a  s. 

iii)  Si  a  >  1 ,  Lf(x)  <  2a  y  Lf-(x)  <  2  -  1/(2a)  en  [a,b],  entonces  la  sucesión  {xn  a} 
es  decreciente  a  s. 

Además,  la  sucesión  {xn  a}  tiene  al  menos  convergencia  cúbica  en  los  casos 
anteriores. 

Demostración  : 

En  primer  lugar  vamos  a  probar,  para  los  tres  casos,  que  la  sucesión  {  xn  a  } 
es  decreciente  y  que  está  acotada  por  s. 

Puesto  que  xn  a-  s  =  Fa(xn_i  a) '  Fa(s),  es  clar0  que  Sl  Fa  es  creciente  en 
[a,b],  aplicando  inducción  y  el  Teorema  del  Valor  Medio  obtenemos  que  la  suce¬ 
sión  está  acotada  por  s.  Entonces  operando  en  la  expresión  (2.8.)  obtenemos 


Lf(x)  2 

F  a<x)  = - 1 - 2  t2  a  (3  -  Lf  -(><))  +  O  -  a)  L  (x)  -3  a]  . 

[2  a  -  Lf  (x)  ]2  f  f 

Cuando  a  <  0  entonces  F’a  (x)  >  0  sí  y  sólo  sí  2a2(3-Lf>(x))  +  (1-  a)  Lf(x)  -  3a 
>  0  y  como  (1-  a)  Lf(x)  >  0,  F’a(x)  >  0  si  2  a^(3  -  Lf-(x))  -3  a  >  0.  Pero  esto  ocurre 
sí  y  sólo  sí  Lf(x)  <  3  -  3/(2a).  Esto  prueba  que  en  el  caso  (i)  la  sucesión  {xn  a} 
está  acotada  por  s.  La  demostración  de  este  hecho  es  análoga  para  (ii).  El  caso 
(iii)  es  similar  teniendo  en  cuenta  que  Lf(x)  <  2a  en  [a,bj. 


f  (x  ) 

'  n- 1 ,  a 


Por  otra  parte  como  x^^-  xn  a  =  f  -  H(a,Lf(xn ^  a))  >  0 

'  n- 1 , cr 


si  y  sólo  si  H  (a,  Lf(xn_i  a))  >  0.  Entonces,  como  H  (a,  Lf(xn_-j  ,a))  -  ^ 
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para  a  <  O  en  R+  y  para  a  >  O  en  (O,  2  a  )  y  además  xn_i  a  >  O,  obtenemos  que 
la  sucesión  {  xn  a}  es  decreciente. 

Por  tanto  {  xn  a}  tiene  un  límite  u.  Tomando  límites  en  (2.8.)  se  tiene  que 
f(u)  =  0,  luego  u  =  s. 

Además,  por  construcción,  todos  estos  procesos  iterativos  tienen  al  menos 
convergencia  cúbica.  # 

A  continuación  nos  planteamos  el  siguiente  problema,  dada  una  ecuación  no 
lineal  (2.1.)  queremos  encontrar  el  valor  de  a  para  el  cual  la  sucesión  obtenida  en 
(2.8.)  converja  a  s  más  rápidamente  .  Para  ello  es  suficiente  estudiar  la  función 


H(  a  ,x )  =  1  +  —  p 

v  '  2a -x  K 


para  un  x  fijo.  Entonces  como 


H(  a  ,x ) 


es  una  función  decreciente  en  a,  obtenemos  el  siguiente  resultado 


Teorema  2.16. 

Sean  ai  ,  i  =  1 , 2,  números  reales,  verificando  ambos  que  ai<  0,  ó  ai  e  (0,1), 
ó  ai  e  [1,  +°°)  y  {xn  a¡),  i  =  1,2  las  sucesiones  decrecientes  a  s  obtenidas  en  las 
condiciones  del  Teorema  2.15.,  con  x  oai=  *oa2=  xo  y  f(xo)  >  0-  Entonces,  si 
al  <  a2  se  sigue  que  xnai  <  xna2  para  todo  n. 

Demostración  : 

Si  al  <  a2  <  0  ,  entonces  H(a1,  Lf(x) )  >  H(a2,  Lfx) ). 

Al  ser  las  funciones  Fa¡,  para  i  =  1,2,  crecientes  probaremos  por  inducción 
clL,e  xn,a1  <  xn,a2- 

Si  0  <  al  <  a2  <  1  ,  como  Lf(x)  <  2a¡  en  [a,b],  para  i  =  1,2,  entonces 
H(a1,Lf(x))  >  H(a2,  Lf(x)  )  y  por  un  argumento  similar  obtenemos  que 
xn,a1  <xn,a2  Para  cada  neN. 

Si  1  <  al  <  a2  ,  la  demostración  es  análoga.  # 

Podemos  deducir  de  este  hecho  que  a  valores  de  a  más  pequeños  la  suce¬ 
sión  {xn  a}  converge  a  s  más  rápidamente. 
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2.3.2.  Notas  prácticas 

Dada  una  función  f,  queremos  elegir  un  valor  de  a  que  nos  dé  la  mejor  suce¬ 
sión  (xn  a},  de  la  familia  (2.8.),  para  resolver  (2.1.).  Para  ello  consideramos  para 
Lf  las  condiciones  que  aparecen  en  el  Teorema  2.15. 

Teniendo  en  cuenta  la  figura  siguiente  podemos  distinguir  las  siguientes 
situaciones  :  M(Lf)  >  3,  M(Lf.)  e  (-°°,  3/2  ),  M(Lf>)e[3/2,  2  )  y  M(Lr)  e  [2,  3]  . 


Fig.29. 

Si  M(Lf)  >  3  existe  solamente  un  a*  e  (-=*>,  0)  tal  que  M(Lf)  =  3  -  3/(2a*)  y 
por  tanto  aplicamos  para  este  a*,  el  proceso  iterativo  (2.8.)  para  resolver  (2.1.), 
siendo  x0  a*e  [a,b]  tal  que  f(x0  a*)  >  0.  Este  proceso  iterativo  converge  a  s  por  el 
Teorema  2.15.  y  es  el  mejor  para  resolver  la  ecuación  (2.1 .)  por  el  Teorema  2.16. 

Si  M(Lf>)  e  (  -oo,  3/2),  existe  solamente  un  a*  e  (0,1)  tal  que  M(Lf>)  = 
3-3/(2a*),  entonces,  como  antes,  aplicamos  el  proceso  iterativo  (2.8. )  para  este  a*. 

Si  M(Lf>)  e  [  3/2,  2),  existe  solamente  un  a*  e  [1 ,  +  °°)  tal  que  M(Lf>)  =  2-1/(2a*), 
y  como  en  los  otros  casos,  aplicamos  el  proceso  iterativo  (2.7.)  para  este  a*. 

En  las  dos  últimas  situaciones,  si  M(Lf)  <  2a*  podemos  elegir  cualquier 
xo,a*  e  [a,t>].  En  otro  caso  ,  si  M(Lf)  >  2a*,  como  Lf(s)=0  existe  c  e  R  tal  que 
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[a,c]c  [a,b]  con  Lf(x)  <  2a*  en  [a,  c].  Y  por  tanto  obtenemos  un  resultado  de  con¬ 
vergencia  local. 

Finalmente,  si  M(Lp)  e  [2,3]  consideramos  el  método  de  Chebyshev  para 
resolver  (2.1 .),  cuya  convergencia  está  asegurada  en  el  Teorema  2.1 . 


Ejemplo  2.3 

Sea  la  función  f(x)  =  x3  +  30x  - 117  ,  entonces  buscamos  el  proceso  iterativo 
más  rápido  de  la  familia  (2.8.)  para  resolver  la  ecuación  no  lineal  f(x)  =  0  en  el 
intervalo  [2,4],  Es  fácil  obtener  que 


Lf(x) 


2  |  60x2  -  702x  -  600 

3  9x4  +  1  80x2  +  900 


Entonces  tenemos  que  Lp  es  una  función  decreciente  en  [2,4],  por  tanto  M( 
Lp)  =  Lp(2)  =  7/4  ,  luego  N  e  [3/2,2).  Entonces  consideramos  M(  Lp)  =  2  -1/(2a*), 
de  donde  obtenemos  a*  =  2. 

Por  otra  parte,  como  60  x3  -  702  x  -  600  <  0  y  9x^  +  1 80  x3  +  900  >  0  en  [2,4], 
se  sigue  que  Lf(x)  <  2/3,  luego  M(  Lf)  <  2/3  y  por  tanto  M(Lf)  <  2a*  =  4.  Así  pode¬ 
mos  tomar  cualquier  x0  2e  [2,4]  tal  que  f(xo2)  >  0,  como  punto  de  salida  de  la 
sucesión  {x0  2},  que  será  la  más  rápida,  que  se  puede  obtener  de  la  familia  (2.8.), 
en  converger  a  la  raíz  de  f(x)  =  0  en  el  intervalo  [2,4]. 
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